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Abstract 

Several substituted coumarins and Schiffs bases have been prepared by the condensation of 
N"(2-ch!oro-5-tri fluoromethyl) phenyl malonamic acid with 2-hydroxy-5-(R)-phenyI azo 
benzaidehyde. 

(Keywords : spectral studies/coumarin derivatives/Schiffs bases/N-(2-chloro-5*tri fluoromethyl) 
phenyl malonamic acid). 


Introduction 

Compounds bearing azo groups and azo coumarins are known to have various biological 
activities^'^. Mishra et al.^ have observed that some coumarin possessed anti-bacterial 
activity against Pseudomonas aeruginosa, E-coli and B-subtilis. 

The present work describes the condensation of N-(2-chloro-5-tri fluoromethyl) phenyl 
malonamic acid with some azo salicylaldehydes to form new coumarins bearing azo group 
at position 6 and new Schiff bases. 

R=RpR2=:R3=H. R=CH3, Ri=R2=R3=H. Ri=CH3, R=Ri=R3=H. R=CH3. R=Ri=R3=H. 

R=C1, Ri=R2=R3=H. R2=C1, Ri=R=R 3=H. R2=Br. R=Ti=R3=H. Ri=R 3=H. RpN02, 

R=:R2=R3=H. 
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Materials and Method 

The general method for the preparation of coumarin and Schiff bases is to heat an 
equimolecular quantity of acid and aldehyde along with a drop or two of pyridine at 105 
toll5“C in an oil bath for 4 h. The mass first melted to a clear liquid which soon set to a 
solid. After completion of heating this was extracted with sodium bicarbonate solution to 
remove unreacted acid and the solid product left was invariably deeply coloured (yellow, 
orange, or mud brown). It was found to be a mixture of correspond! ng coumarins and Schiff 
bases. The Schiff base was easily separated from the coumarin being soluble in ethanol, 
while the coumarins remained insoluble. The coumarins were purified be repeated washing 
with boiling solvents such as ethanol, methanol and acetone. The Schiff bases were 
recrystallised from ethanol. 

The melting points were determined using sulphuric acid bath. 1. R. Spectra in Nujol 
were recorded on a Perkin Elmer spectrophotometer. 'H-NMR spectrum were recorded on 
Brucker 300 MHz and TMS as an internal standard (chemical shifts in 5 ppm). 

Results and Discussion 

The molecular formula, M.P., % yield, percentage of nitrogen, ‘ H-NMR and HRMS 
of various coumarins and Schiffs bases are reported in Tables (1 and 2). 

The filter paper disc method"* was used for screening the compounds for antimicrobial 
activity. Standard size discs of Whatman filter paper No. 3 with a diameter of 6.5 mm were 



Table 1 Coumarins obtained by the condensation of N-(2-chloro-5-trifluoromethyI) phenyl malonamic acid with 2-hydroxy-5-(R)-phenyl i 
benzaldehyde. 
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sterilised by dry heat at 140“C for 1 h and were saturated with the solution of compound. 
These were air dried at room temperature to remove any residual solvent, which interfere 
with the determination. The discs were then placed on the surface of the sterilised solidified 
culture medium that had been innoculated with the test organism (using a sterile cotton 
swab) and air dried to remove the surface moisture. The thickness of the culture medium 
was kept equal in all petridish (20 mi). For each of the experiment a petridish innoculated 
withthe organism but without disc containing the compound was taken and incubated as 
control. These were incubated at 28 ± VC after which the zone of inhibition or depressed 
growth were measured. 

The data of antibacterial screening is reported in Table 3. The effectiveness of these 
compounds were tested against S. aureus, E. coli, Pseudomonas and Klebsiella. 


Table 3 - Anti-bacterial activity screening of coumarin. 


Name of the coumarin and 



Zone of Inhibition (mm) 


dilution 

S. aureus 

E. coli 

Pseudomonas 

Klbsiella 

6-(4-chloro)pheny! 

azo-coumarin-3-carboxy- 

(2-chIoro-5-tri 

fluoromethyl) 

phenyl amide 





1 :5 

R 

R 

R 

R 

1 : 10 

R 

R 

24 

R 

1 :20 

R 

R 

R 

R 

1 :50 

R 

18 

22 

R 

1 : 100 

R 

R 

R 

R 


R = Resistant 

None of the other coumarins showed antibacterial activity. 

None of the schiff bases showed antibacterial activity. 


The purity of the compounds were checked by T.L.C. 5% benzene/methanol for Schiff 
bases and 10% benzene/methanol for coumarins. 

IR spectrum of 6-(2, 3-dimethyi) phenyl azo-coumarin-3-carboxy“(2- chioro-5-tri 
fluoromethyl) phenyl amide shows absorption at 3400 cm”^ (-NH stretching vibrations), 
1720 cm~^ (C=0 of secondary amide), 1660 cm"^ (C=C stretching), 1580 cm”^ (N=N 
stretching), 1320 cm”^ (-CH3 bending), 1110 cm""^ (aromatic character), 1070 cm”^ (C-F 
linkage), 810 cm'^ (C-Cl linkage). 
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Materials and Method 

SMZ (m.p. i 69 *C) used in the present study was procured from Matta Remedies, 
Ludhiana (India). All the metal salts i.e. AgN03, CdCl2, HgCl2, C0CI2.6H2O, VOSO4.H2O 
and U02(N03)2.6H20 employed were reagent grade chemicals. All the complexes were 
prepared as follows^, A solution containing 5 mmol of SMZ in aqueous sodium hydroxide 
(pH 8-iO) was added to a solution of 2,5 mmol metal ion namely VO(II), U02(II), Co(II), 
Hglll), Cd(n) and Ag(I). The precipitates that formed immediately or after a few minutes 
of stirring were filtered, washed with distilled water and dried in a vacuum desiccator. 

Method of Analysis : Melting points were determined with a Thiele’s apparatus using 
open capillar>^ tube. Carbon, hydrogen and nitrogen were analysed with a Perkin Elmer 
240 C elemental analyser. Sulphur and metals were estimated by the standard methods 
reported in the literature^'^. Infrared spectra of the ligand and its metal complexes were 
recorded in nujol mull on a Perkin Elmer 1430 double-beam spectrophotometer. 

The magnetic susceptibilities of the metal complexes prepared were measured with 
Gouy’s method at room temperature using Hg[Co(NCS)4] as calibrant. Conductivity 
measurements w^ere carried out on bridge type CM 82 T conductometer. Differential scanning 
calorimetric (DSC) studies were carried out on a Stanton Rederoft differential scanning 
calorimeter using AI2O3 as a standard and zinc as reference. The heat of fusion of the metal 
complexes was determined using well known standard formula^. 

H} KLAJm,)M, 

where = heat of fusion of Zn, 7.28 kJ mor*’ = heat of fusion of sample, K = 
calorimeter constant, mz„ = mass of zinc, = mass of sample, = area of the Zn peak, 

= area of the sample peak, M 2 „ = atomic mass of Zn, = molecular mass of sample. 
The area of each DSC peak was determined by paper weight method. 

Results and Discussion 

The analytical and nonspectral physical data of the ligand and metal complexes have 
been listed in Table 1. The stoichiometry (M ; L) of the complexes as deduced from 
elemental analysis has been found to be 1:2 in case of Cd(n), Co(II), VO(II), UOjCII) and 
Hg(ll) complexes while 1:1 in Ag(I) complex. The stoichiometry of the complexes also 
conforms to the thermogravimetric and differential thermal analysis data". These complexes 
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Fig. 1 - Structure of sulphamethoxazole. 


have high melting point and are insoluble in water and common organic solvents namely 
chloroform, ether, benzene, carbon tetrachloride and ethanol. The molar conductivity value 
of complexes in DMF of Ag(I), 19.3; Cd(II), 32.0; Co(II), 45.5 and Hg(II), 4.5 ohm ^ cm^ 
mo!”^ show these compounds to be non-electrolytes. 

IR spectra : The sulphamethoxazole moiety is expected to behave as chelating ligand 
with the metal ions. Just like its other analogues, the chelation^^’^^ occur mainly through 
(i) nitrogen of NH 2 and sulphonyl oxygen of - SO 2 NH groups, (ii) nitrogen of -NH 2 as 
well as --SO 2 NH groups, (iii) sulphonyl oxygen and nitrogen atom of -SO 2 NH group and 
(iv) sulphonamide nitrogen and ring nitrogen of isoxazole ring. The NH 2 stretching, NH 2 
bending, NH and SO 2 stretching (- SO 2 NH), C-N stretching (- SO 2 NH and NH 2 ) and C=C 
stretching (phenyl ring), 0-H stretching (water) vibrations of the ligand and its metal 
complexes has been assigned^’^"^. and recorded in Table 2. A perusal of IR data in Table 2 
reveal that -NH 2 and - SO 2 NH groups of SMZ are involved in coordination since the 
region corresponding to the vibrational bands of these groups are affected on complexation. 
The bands of medium intensity of SMZ at 3480 and 3387 cm“^ due to asymmetric and 
symmetric -NH 2 stretching vibrations of the amino group except Co(II) and Cd(II) are not 
influenced significantly on complexation. This suggests the participation of nitrogen atom 
of “-NH 2 group in Co(II) and Cd(II) complexes. This is also in accordance with the 
observations of Narang and Gupta^^ on the metal complexes of the related sulpha drugs 
i.e. sulphanilamide, sulphaguanidine, sulphathiazole, sulphamerazine, sulphadiazine and 
sulphapyridine. The nitrogen of the isoxazole ring of SMZ is not involved in coordination 
in the present complexes since the G=N stretching frequency is not affected at all on 
complexation. Further in the Co(II) and Cd(II) complexes, the significant shift in the SO 2 
symmetric stretching frequency and with change in intensity also infer the involvement of 
oxygen of sulphonamide group (Table 2). This conforms to the observations of Kanagaraj 
and Rao^ on Cr(III), Mn(II), Ni(II) and Cu(II) complexes. In view of it, we would like to 
point out that the observations of Kanagaraj and Rao^ on the basis of infrared and NMR 
spectral data that absence of (i) vibrational band at 3300 cm"^ (due to sulphonamide N-H 
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medium; w - weak; vw - very weak; sh - shoulder. 
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stretching of the free ligand) in the metal complexes and (ii) N-H proton signal at 
11 . 195 , ppm (which disappears on deuterium exchange) in the zinc complex due to 
deprotonation seem to be erroneous because the complexes were prepared in the alkaline 
medium (pH 8-10) and the pKa of the ligand is 6^^. In this pH range, the ligand itself 
undergoes deprotonation. Hence the absence of N-H band at 3300 cm~^ in the infrared and 
proton signal at 11.195, ppm in the NMR spectra cannot be fixed as criterion for the 
deprotonation of sulphonamide proton of SMZ on complexation and hence the binding of 
metal ion Co nitrogen of the SO 2 NH group as proposed by the authors. Most probably 
inadvertently Kanagaraj and Rao^ has misquoted Gogorishvili et In sulphonamide 
ligand, the latter have suggested the binding of metal ions to the ligand through the nitrogen 
of the isoxazoie ring since the C=N band of the ligand at 1600 cm"^ appears as doublet 
(1598-1605 and 1587-1595 cm“^) in the metal ion complexes. At no stage, they have 
mentioned that the missing of the band at 1380 cm"^ and a positive shift of the band at 
1460 to 1480 cm \ indicate that the involvement of isoxazoie ring in the coordination of 
the ligand (possibly through nitrogen) to metal ion. In fact, they have assigned a band at 
1380 cm'^ to S=0 asymmetric stretch in the IR spectra of sulphanilamide. This band appears 
in the range 1378-1385 cm ^ in complexes. Thus, the isoxazoie ring vibrations of SMZ at 
1460 and 1380 appearing at 1490 and 1480 cm"-^ in its copper complex may not be 
fixed as criterion as suggested by Kanagaraj and Rao"^ for the coordination of the metal ion 
to the ring nitrogen since the C=N stretching vibration of the same very ring remains 
unaltered in the complexes. 


Just like Co(II) and Cd(II) complexes of SMZ, Ag(I), VO(II), U02(II) and Hg(in also 
bind to SMZ through oxygen of the SO 2 NH group since S=0 (symmetric/asymmetric or 
both) are influenced on going from ligand to the respective complexes (Table 2). However 
in the case of UO^fll) and Hg(II) complexes, nitrogen of SO,NH is also involved in 
coordination as is evident from the C-N stretching frequency shift and intensity changes 
in passing from the ligand to the complexes (Table 2). The shift is found to be of the orLr 
of 4- 2 and 17-19 cm’ in the Hg(II) and UO^dl) complexes respectively. All the metal 
complexes except that of Co(II) were found to be diamagnetic. Though there are various 
factors 1 e stenc, pH, electronic state of the metal, which govern the extent of co-ordination 
with various metal ions, no single factor in this case can explain 
the extent of co-ordmation of SMZ with different metal ions because of their diverse naC 

SMz'^nTrl^f ^ characteristic parameters of the DSC curves of 

SMZ and its metal complexes are recorded in Table 3. Heat changes associnter) u,i,h „ 
steps of dissociation of the metal complexes as deduced from efararr' ' 

DS^urvK are also recorded in Table 3. The bread oharecer of some of the DSC Tab 
as obsereed .. vanoas complexes sogpests tha, a number of reaefions JoceZg 
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Table 3 - Characteristic parameters of DSC curves and heat of fusion of SlRlZ and its metal complexes. 


Compound 

Temperature (“C) 

Peak 



A5 


range 

temp. (“O 

(mg) 

(kJ mor') 

(Jmor'K"') 

C 10 H 11 N 3 O 3 S 

17(f - 198 (s, endo) 

180 

2.20 

38.18' 

68.18 


230-290 (s, exo) 
325-610 (b, exo) 

265 

8.80 

152.73 


[Ag(SMZ)H20] 

249"-310 (s, exo) 
500-570 (b, exo) 
570-620 (b, exo) 
650-720 (b, exo) 

280 

9.00 

233.22 

446.78 

[Cd(SMZ)2(H20)2] 

160-20Cf (s, endo) 

189 

3.00 

134.23' 

280.20 


280-309 (s, exo) 

295 

0.35 

15.66 



310-360 (s, exo) 
500-670 (b, exo) 

330 

2.70 

120.80 


[V0(SMZ)2(H20)2] 

100-130 (s, endo) 

110 

1.82 

34.13 

133.79 


175^-205 (s, endo) 

190 

1.44 

59.94'= 



240-285 (s, exo) 
360-525 (b, exo) 
605-700 (b, exo) 

265 

7.68 

319.68 


[U02(SMZ)2].H20 

100-135 (s, endo) 

120 

0.70 

38.11 

157.28 


177^-210(8, endo) 

190 

1.30 

70.78^= 



235-350 (s, exo) 
440-700 (b, exo) 

290 

6.00 

326.68 


[Co(SMZ)2(H20)2].H20 

152^-215(8, endo) 

185 

3.95 

167.13" 

393.24 


245-305 (8, exo) 

' 500-570 (b, exo) 
660-700 (b, exo) 

279 

5.00 

211.56 


[Hg(SMZ)2(H20)2] 

198-230“ (s, endo) 

210 

0.40 

20.31" 

40.37 


250-315 (8, exo) 

280 

4.30 

218.38 



600-700 (b, exo) 

s-symmetncal: b-broad; exo-exothermic, endo-endotheimic, Hrheat of fusion; AG=0; AA*’-area 
under thp DSC curve; AA„- 1.3 mg; *iiielting temperature. ■ 

simultaneously with the formation of unstable intermediate species in the decomposition 
process^ as expected. The value of (with superscript c) in column five (Table 3) 
corresponds to the heat of fusion of the respective systems while the second value 
corresponds to some sort of decomposition of the respective systems. The heat of fusion 
reflects on the thermal stability of a metal complex in a particular state . Thus the 



202 


B.S. SEKHON et al 


trend of the thermal stability of the metal complex on the basis of (Table 3) is : Co(II) 
> Cd(ll) > UO.di) > VO(II) > Hg(II). 


Conclusions 

On the basis of spectroscopic and differential scanning calorimetric studies on the 
prepared metal complexes of SMZ, the following aspects of the behaviour of SMZ are 
demonstrated. SMZ acts as (a) monodentate ligand binding the Ag(I) and VO(II) through 
the oxygen of the SO 2 NH group, (b) bidentate ligand binding the metal through nitrogen 
of NH 2 and oxygen of SO 2 NH group in case of Co(II) and Cd(II) complexes while via 
both oxygen and nitrogen of SO 2 NH group in 1102(11) and Hg(n) complexes, (c) the absence 
of N-H band at 3300 cm”^ in the infrared and proton signal at 1\195, ppm in the NMR 
spectra cannot be fixed as criterion for the deprotonation of sulphonamide proton of SMZ 
on complexation and hence the binding of metal ion to nitrogen of SO 2 NH group as proposed 
by Kanagaraj and Rao^, (d) the thermal stability of the metal complexes follows the order 
; Co(II) > Cd(ll) > U02(n) > ¥0(11) > HgdI). 
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Abstract 

The title tercopolymer (S-HQOF) have been prepared by refluxing a mixture of 8-hydroxy quinoline 
(8-HQ), oxamide (O) and formaldehyde (F) in presence of 2M HCl as a catalyst for 4-5h at 127 (± 

1 “O with varying molar ratio of reactants. The fourdiffemt tercopolymer thus obtained using different 
molar ratios are ail amorphous in nature with four different-shades of light brown colour. They have 
high melting points and exhibit high thermal stability. The tercopolymers were analysed for C, H and 
N. Their purity was established by thin layer chromatography (TLC). The viscosity measurements, 
in dimethylformamide (D.M.F.) have been carried out with a view to ascertain characteristic funcdon 
and constants. The molecular weight of tercopolymers were determined by non-aqueous 
conductometric titrations against alcoholic KOH in DMF and also ultra-violet, infra-red and nuclear 
magnetic resonance spectral studies have been undertaken for synthesized resins to elucidate their 
structure. 

(Keywords : 8 HQOF/tercopolymer) 


Introduction 

A large number of synthetic resins derived from hydroxy and amino compounds have 
attracted the attention of many research workers due to their versatile use as ion 
exchangers, photographic binders^ and thermal stabilizers^'^ etc. 

Chelating resins containing oxine group have been synthesized either from oxine, 
formaldehyde and resorcinol by condensation^’^ or by diazotization of poly (amino styrene) 
resins followed by coupling with oxine^. Their properties were later improved by Pennington 
and Williams, and also by Parrish and Stevenson^^. Aristov and Kostaintinov reported 
preparation of chelating ion-exchanger by copolycondensation of 8-hydro xyquinoline or 


*Paper presented at “The National Academy of Sciences, India, 68th Annual Session, CDRI, Lucknow. 

(Oct. 23-25, 1998) 
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phenol derivatives like amino phenol, P-resorcyclic acid or resorcinol with formaldehyde^ \ 
Perohit, ei al prepared a resin from 8-hydroxyquinoline, hydroquinone and formaldehyde^^. 
8-Hydroxy quinoline, resorcinol, formaldehyde and 8-hydroxyquinoline, resorcinol, 
salicyciic acid, formaldehyde, chelating resins were prepared and studied for its 
ion-exchange properties by Kim Don et al}^ Kapadia and co-workers have prepared many 
resins containing 8-hydroxyquinoline by using 8-hydroxyquinoline and poly (vinyl acetate) 
in glacial acetic acid^^ 8-hydroxyquinoline and poly (vinyl alcohol) in DMF^^, 
8-hydroxyquinoline and epichlorohydrine with dimethyltetramine^^, and 
8-hydroxyquinoline and acetaldehyde with thiourea^^. 

However, literature sur\^ey revealed that no tercopolymers have been synthesized from 
8-hydroxyquinoIine, oxamide and formaldehyde. Therefore, we have undertaken the 
synthesis and characterization of these tercopolymers and the results of our studies are 
reported in this paper. 


Materials and Method 

All chemicals used were of AR Grade. The solvents like dimethyl formamide and 
dimethylsulphoxide were double distilled under reduced pressure. 

Synthesis ofS-hydroxyquinoline-oxamide-fomialdehyde (1) tercopolymer : A mixture 
of 0.1 mole 8-hydroxyquinoline, 0.1 mole oxamide, and 0.2 mole formaldehyde in 100 ml 
of 2M HCl was refluxed on oil bath, at 127 (± TC) for 4-5 h with occasional shaking. The 
contents were then cooled to room temperature and neutralized with dil. NaOH. The greenish 
yellow resinous product so obtained was filtered and repeatedly washed with cold water, 
hot water and alcohol and dried in air. The powder was then purified by dissolving it in 
8% NaOH solution, filtered and reprecipitated by gradual dropwise addition of ice cold 
1 : 1 (v/v) concentrated HCl and distilled water with constant and rapid stirring to avoid 
lump formation. The tercopolymer resin so obtained finally was filtered, washed several 
times with hot water, dried in air, powdered and kept in vacuum desiccator over silica. 

Similarly, other tercopolymer resins viz. 8-HQOF (11), 8-HQOF (III) and 8-HQOF (IV) 
were synthesized by varying the molar proportion of the starting materials i.e. 
8-hydroxyquinoline, oxamide and formaldehyde in the ratios 2:1:3, 3:1:4 and 4:1:5 
respectively. 

Characterization i 

Methods : The purity of each tercopolymer sample was checked and confirmed with 
thin layer chromatography. Elemental analysis of all the tercopolymers for carbon, hydrogen 
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and nitrogen were carried out in CARLO ERBA STUM, DP 200 at CDRI Lucknow. 
Molecular weights of the resins were determined by non-aqueous conductometric titration 
in DMF. The viscosity measurements were carried out in DMF at 35°C using Tuan-Fuoss 
viscometer. Electonic absorption spectra of all the newly synthesized tercopolymer resins 
were recorded in DMSO on Shimadzu double beam spectophotometer. The Infrared spectra 
of all the newly synthesized tercopolymer resins have been scanned on Perkin-Elmer 983 
spectrophotometer in KBr pellets in the wave number region of 4000-400 cm“^ Proton 
NMR spectra of all the newly synthesized tercopolymer resins have been scanned on 
Varian-EM 360 A 60 MHz proton, NMR spectrophotometer using deuterated 
dimethylsulphoxide (DMSO-D^) as a sovent. 

Results and Discussion 

All the tercopolymers were light brown in colour and were completely soluble in DMF, 
DMSO, HCl and aq. NaOH. The details of synthesis and analytical data of all the 
tercopolymers are presented in Table 1. 

Molecular weight determination (Mn) : The results of conductometric titrations were 
used to determine number average molecular weight (Mn) of the tercopolymer samples 
using alcoholic KOH as titrant^^’^^. From the plots of specific conductance against 
milliequivalents of titrant (ale. KOH ) added, the first breaks and the last breaks were noted. 
The degree of polymerization (DP) of the tercopolymer samples was obtained from the 
values of the breaks. Table 2 shows that the molecular weights increase as proportion of 
8-hydroxyquinoline increases in the tecopolymer resins^^. 

Viscosity measurements : The measurements of viscosity have been carried out in DMF 
at 35^*0 using Tuan-Fuoss viscometer. The reduced viscosity versus various concentrations 
is plotted for each set of data. The intrinsic viscosity [q] has been determined from the 
following equations. 

Huggin’s equation 

%p/c = + *1 c ( 1 ) 

Kraemer equation 

In (Ti/q = [r|] - [r\f C (2) 

The intrinsic viscosity of 8-hydroxyquinoline-oxamide-formaldehyde tercopolymer 
determined from both the plots using eqn. (1) and (2) is identical and it increases with the 
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Fig. 1 ~ Structure of 8-HQOF - 1 tcrcopolymer 



Fig. 2 - Structure of 8-HQOF - II tcrcopolymer 



Fig. 3 Structure of 8— HQOF — III tercopolyitier 


r 



Fig. 4 - Structure of 8-HQOF - IV tcrcopolymer 



SYNTHESIS & CHARACTERIZATION OF TERCOPOLYMERS (8-HQOF) 


209 


Table 1 - Synthesis and analytical data of 8-hydroxyquinoline - oxamide ~ formaldehyde (8-HQOF) 
tercopolymers. 


Tercopolymer 


Reactants 


Colour 

Yield 

% 

C % 

Obs 

(Cal) 

H % 

Obs 

(Cal) 

N% 

Obs 


8-HQ 

mol 

Oxamide 

mol 

HCHO 

mol 



(Cal) 

8-HQOF - I 

1 

1 

2 

Light 

brown 

65 

61.10 

(61.88) 

4.30 

(4.34) 

16.34 

(16.42) 

8 - HQOF - II 

2 

1 

3 

Light 

brown 

67 

67.20 

(67.78) 

4.34 

(4.81) 

13.32 

(13.62) 

8 - HQOF - III 

3 

1 

4 

Light 

brown 

72 

69.30 

(69.66) 

4.38 

(4.71) 

12.30 

(12.64) 

8 - HQOF - IV 

4 

I 

5 

Light 

brown 

80 

70.88 

(70.52) 

4.40 

(4.52) 

11.54 

(11.86) 


increase in molecular weights of tercopolymer samples under investigation (Table 2). 
Similar trends have been reported by earlier workers^^’^^ 

UV-visible spectra : UV-visible spectra of all the four 8>hydroxyquinoline-oxamide- 
formaldehyde tercopolymers show two characteristic bands at 330 nm and 270 nm. The 
intense band with some dimnished fine structues at 330 nm is due to (n n*) allowed 
transition of the benzenoid system of heteroaromatic compounds (quinoline). The 
bathochromic shift from the basic value of 308 nm may be due to presence of phenolic 
hydroxy group (auxochrome). The second less intense band at 270 nm is due to (n tx*) 
forbidden transition in saturated aliphatic carbonyl compounds . The appearance of 
phenolic hydroxyl group as well affects the rise of € value, i.e. hyperchromic effect. 
The log value gradually increases in the order : 8-HQOF (I) (log 0.41) < 
8-HQOF (II) (log 0.45) < 8-HQOF (III) (log 0.46) < 8-HQOF (IV) 
(log 0.55) which is due to introduction of more and more chroraophores (hetero 
aromatic system) and auxochromes (phenolic OH group) in the repeat unit of the 
tercopolymer resins^^. 


IR-Spectra : IR-Spectra of all the tercopolymer resins recorded in KBr matrix resembled 
each other and assigned to the tercopolymer on the basis of nature and reactive positions 
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Table 2 - Molecular weight determination, viscometric data and IR frequencies of 
8-hydroxyquinoline-oxamide-formaldehyde (8-HQOF) tercopolymers. 

Tercopolymers Number average Intrinsic viscosity Important IR frequencies 

8-HQOF mol. weight (Mn) [T]j (dl/g) — ___ 





Wave number 

Assignment’s 

S-HQOF-I 

5140 

8.50 

3371 (mbr) 

- OH Stretching 




2930 (m) 

- CH 2 Stretching 




3050 (mbr) 

> NH. Stretching 




1621 (s) 

>C=0 Stretching 

8 -HQOF~n 

5970 

10.00 

1418 (m) 

-€=N- Stretching 




1230 (m) 

—OH in plane bending 




1579 (m) 

Aromatic ring 

8 -HQOF-III 

10,849 

11.50 

2850 (w) 

-CH 2 Stretching 




1466 (m) 

“CH 2 Scissoring 




1267 (m) 

5,7,8 - trisubstituted 
aromatic ring coxinet 

8 -HQOF-IV 

11,070 

13.00 

70 (m) 

-CH 2 rocking 

m - medium; s - strong; w - 

weak; br - broad. 



^ S-hydroxyquinoIine-oxamide-formaldehyde ( 8 -HQOF) tercopolymers (in 

Nature of proton assigned 


Observed chemical shifts. (5) ppm 

Expected (5) 


8 -HQOF-I 

8 -HQOF-n 

ppm 

8 -HQOF-IIl 8 -HQOF-IV 

Proton of phenolic - OH 

8.75-8.65 

8.4 - 8.65 

8.85 - 8.65 8.85 

-8.65 8.0- 12.0 

involved in hydrogen bonding 


Aromatic proton (Ar - H) 

8.33-8.14 

8.49 - 8.23 

8.36 - 8.24 8.44 

- 8.24 7.3 - 8.8 

Proton of “ NH - bridges 

7.34-6.65 

7.49 - 7.08 

7.51-6.14 7.53- 

- 6.11 5.0 - 8.5 

Amido proton Ar-CHaNH - 
CO ~ linkage 

4.10-3.39 

4.64 - 3.88 

4.53 - 3.36 4.67 - 

- 3.63 3.5 - 6.00 

Methylene proton of Ar - CH 2 
“ Ar - linkage 

- 

3.48 - 2.49 

2.81-2.48 2.50- 

- 2.48 2.0 - 3.0 
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of the monomers. A broad peak in the range 3360-3377 cm~^ was assigned to hydrogen 
bonded v A broad band observed around 3000 cm“^ may be assigned to v NH of 
the oxamide monomer^"^’^^. A sharp and strong band observed at around 1621 cm“^ may 
be assigned to v C=0 of the tercopolymer“^’^^. A sharp peak around 1500-1580 cm“^ may 
be assigned to aromatic skeletal ring breathing modes^^. The inflections at around 2800, 
1460, 1260 and 700 cm~^ due to stretching, scissoring, wagging and rocking modes of 
-CH 2 -- suggested the presence of methylene bridges in tercopolymer^^’^^. 

NMR Spectra : The three 8-hydroxyquinoline - oxamide - formaldehyde 
tercopolymers-II, III and IV exhibit signals in the region 2.5 to 2,8 5 (ppm) which are 
attributed to methylene proton of Ar-CH 2 “-Ar linkage*^. A broad signal at 3 to 45 (ppm) 
may be assigned to methylenic proton of Ar~CH 2 -NH^O linkage^^. The presence of broad 
signal may be attributed to the presence of -NH-bridges^^ The signal around 8.1 to 8.3 
and 8.6 to 8.95 (ppm) are attributed to proton of aromatic ring (Ar-H) and phenolic-OH 
involved in hydrogen bonding respectively^^’^^. 

On the basis of the nature and reactive positions of the monomers, elemental analysis, 
IR and NMR spectra and molecular weight, the following structures (Fig. 1^) have been 
suggested for the tercopolymers. 
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Abstract 

The approximate model introduced by Leung in 1984 was rediscussed by Gupta and Srinivasan, who 
point out that the model refers to a uniform distribution of B-trees, i.e., a distribution in which all 
possible configuration of B-trees with the same number of keys are equally likely. To be more precise, 
if /is the minimum storage utilization of each node, then each node can contain at least Idfkeys and, 
obviously, at most 2d keys; Gupta and Srinivasan observe that the problem of finding the number of 
different B-tree configurations having an order d, n nodes and r keys is equivalent to the restricted 
occupancy problem when n different cells can contain more than v of F like objects, where P = 
r - 2dnf and v = {l-f)2d. 

Unfortunately, this statement is not precise and Gupta and Srinivasan’s model can not be considered 
to be a true model for B-trees. It makes sense to reconsider Gupta and Srinivasan’s model from this 
point of view and in fact we obtain more precise results on storage utilization. By the way, from a 
mathematical point of view, these results conclusively show that the model cannot be considered a 
realistic B-tree model. 

In the present paper, we want to show that it is possible to find the generating function for the sums 
Np (n, v) and ^/z nNp {n, v). 

This allows us to give either approximate or asymptotic estimates for the expected value E{u) and 
thus solve the storage utilization problem for this model. To this purpose we use the method of Riordan 
arrays. Riordan arrays constitute a flexible algebraic approach to the solution of combinatorial sums, 
to find either their closed form or their asymptotic estimate, and seem to have many applications to 
combinatorial analysis and to analysis of algorithms. 

(Keywords : B-trees/Riordan array) 


Introduction 

5-trees are one of the most frequently used structures for storing information in a two 
level storage. Since there introduction in the early 1970’s, they have been widely studied 
together with their main memory companion, the 2-3 trees, and therefore a vast literature 
exists, in which their performance is analyzed. Just to give some examples, we quote here 
the papers by Beaza-Yates^ Chu and Knott^, Mishra^'"^, Quitzow and Klopproge^ and 
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Wright^ which all refer to a distribution in which key permutations are all equally likely. 
This problem is quite hard, and therefore several approximate models have been proposed, 
in order to simplify the analysis, especially concerning their “occupancy problem”, i.e, 
how much space is necessary to store a given number of keys. If we have r-keys, stored 
in a « nodes, each containing a maximum of 2d keys (d is the order of the tree), then the 
ratio r/(2dn) is called the load factor for the B-tree and is the quantity of main interest. 

Preliminaries 

Gupta and Srinivasan’s’ paper shows that, if Np (n, v) represents the number of B-trees 
with p keys distributed over/i nodes such that each node can have a maximum of v vacancies, 
then 


(«. v) 


n 


S (- 1/ 


f"' 


^ P-kv-k+n-i ' 

k 


P-kv-k 

\ J 

K y 


( 1 ) 


In this model, the number n of nodes varies between = r/{2d) and = rl{2df) and 
the expected storage utilization is : 


r 5] Np (n, v) 

£(„) = —fps!! = rl2dE{\ln). 

2d 2) nNp{n,v) 


This definition of E{\ln) is slightly different from Gupta and Sriniva.san’s’ but it is more 
appropriate from a probabilistic point of view. The sums can be extended from 0 to » 
because by (1) the value of N («, v) is zero when n is outside the ranee n n 

H tD rnin ***• 


Riordan Arrays 

The technique we use to solve the occupancy problem relies on the concept of Riordan 
array . A Riordan array \p, ksN, k<pj is an infinite, low triangular array defined by 
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Table I - E{u) = approximate storage utilization. 


D 

5 

/= 1/2 

E{u} 

/= 

i 

2/3 

E{u) 

it 

3/4 

E(u) 

6 

90.796544 

67.72% 

0.849077 

80.64% 

0.865395 

86.12% 

12 

0.859757 

66.71% 

0.894105 

80.28% 

0.910745 

85.95% 

24 

0.908698 

65.50% 

0.930175 

79.80% 

0.940164 

S5J2% 

48 

0.943159 

64.21% 

0.956165 

79.25% 

0.961960 

85.43% 

96 

0.965760 

62.92% 

0.973493 

78.64% 

0.976783 

85.10% 

192 

0.979948 

61.70% 

0.984418 

78.01% 

0.986260 

84.76% 


means of two analytic functions h{t)), in the sense that the generic element < 3 ^^ is 
given by 


ca^ j, = [f] bit) (t(k(t))^) 


[tp] being the “coefficient of’ operator. 


In other words, b{t) {t(h)(t))k, k = 0.1 , is the generating function of the kxh column of 

the Riordan array. Riordan arrays allow us to solve combinatorial sums like ^ k 8k 

k 


terms of generating functions according to the following rule : 


G 


X ^p^kSk 

k=o 


= bit) it (hit)) 


where the symbol G denotes the generating function operator i.e. 
git) = G {g^licg N 


( 2 ) 


= G {g,} 



it=0 
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is the generating function of the sequence 

In particular, the arrays row sums and weighted row sums in which 
g" = k, have the following generating functions : 


G 



..Ht) 


G 


X kap,k]^ 


t b(t) h{t) 


we wish to point out that in formuJa (1), 


^pjc 


(p+/ 2 -l~/:(v+l) 
p-k (v+1) 


is the Riordan array D = (l/(l-r)”, t') and by a simple application of rule (2) with 
= i\-t)n. The generating function of binomial coefficients with alternating signs, we obtJn' 


= X 


^ {p-kn-l-k{v+\)' 


["] 

p-k{v+\) j 


1 

/J 


(-]/. 




(l-r)" > 


= [f] 


l-r' 


v+r 


(i-O 


= [r^^] ( 1 + /+...+ 

+ ;2{l-y)</y. 
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This proves that, when 2d/and 2i\-f)d are integer numbers, then (n, v) is the genetic 
element of the Riordan array* defined by b(t) = 1 and h{t) = (i +?+.,. 

This immediately gives 

S = 2 


- (i-r2#(i + r + .... + r2(i-/)^)’ 

ly = n Np {n, v) 


The most important cases are 

(i) /= 1/2 

This corresponds to (B)-trees whose analogue is the standard definition of R-trees. 

(ii) /=2/3. 

This corresponds to (B)-trees with left [right] rotations. In the B-tree analogue, when a 
node is full, its keys are “rotated” into the left [right] sibling. Node splitting is performed 
only when the two nodes are full. 

(iii) /=3/4. 

This corresponds to (B)-trees with left and right rotations. In the B-tree analogue, when a 
node is full, its keys are “rotated” into the left or right sibling, according to whichever 
has room. Node splitting is performed only when the three nodes are full. 


(iv)/= 1/3. 
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Thib also corresponds to {B)“trees, since the circles of radii exist between 1/2 and 2 and 2 
in B-trees. In any case, we can always, choose d in such a way that both 2d/and (21-/) d 
are integers. The expected storage utilization is given by 


E{ii) = 


rs 

2dw * 


If we now set 


P(t) = 1 (1 




= 1 ~r 


2df 


i_,2^(l-/)^l 


1 


Q{!) = 1 + ... 


It is obvious that the asymptotic values of S and W are connected to the zeroes of P{t) of 
minimal modules for a detailed explanation of singularity analysis for generating functions. 
By applying the results about the location of the zeros of polynomials, we find that the 
zeros of P{t) lie in the annulus bounded by the two concentric circles of radii 1/2 and 2, 
centered at the origin. We observe that, by adding the root r = 1 to P(t), we get the equivalent 
equation. 


= 0 (3) 

This implies that only two real roots exist for P(t) = 0, and J, with 1/2 < 5 < 1 and -2 
< 5 < -1 and multiplicity 1. 

The other 2d-2 roots are all complex. Since the coefficients of the Taylor expansion 
of P(t) ^ around r = 0 are all positive (i.e. our {n, v), should be the root of minimal 
modules) and thus we have : 

1 A A 
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ew . B ^_c 


P{tf {l-i/sf l-t/s 




nr-2df 


r-2df 


for suitable constants A, B and C which are easily computed in terms of S 


A = --J— , 5 = -^ 

•fp‘ is) ' (sp‘ (5))^ 

sp{sf 

A’s and B's values are the most important ones, and can be described in further detail by 
observing that 


Q(s) = \+s + s^+...+s^^'--f^ 


1 _ ^M1-A+1 
( 1 -^) 


^ _ s^f. (1 - - 25(1 ( 1 - 5) 

( 1 - 5 ) ( 1 - 5 )^ 


-2dfs^^~^ .(1 
l-s 


l-s 


1 ' Qjs) 

55 ’ ( 5 ) ’ ( 55 * ( 5 ))^ 


A = - 
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r = (^) - Q' is) P' (s) 

spisf 

Here we use the identity ^ ^ 

because s is a root of P{t). We can now compute the expected values of storage utilization 
as follows ; 


Eiu) ~ 

s-‘‘^2d{B{r-2df+ 1) + C) 

^ ^ ^ -P'(s) 

~ 2ds^^B 2ds-‘‘^'^ Q(s) 

We set r — > CO, i.e. we assumed that our (B) trees contain a large number of keys. By 
inserting the values of P'(s) and Q(s) in this expression, we have 


^ 2d(l-s) 

which is our most significant result. It gives us the expected value of storage utilization in 
terms of the unique real solution of P(t) = 0 sJ. 112 < s < 1. This result gives us some 
further information : by construction, we know that E{u) > /, and this means that 
1 - {2il ( I -f) + 1 + s^d > 0. This provides us with an upper limit for .s : 


" - did/ {2d (I -/} + !} 


If we set e'‘ = 2d{l-f)+ I or 


h = In i2dil-f) + l) 


In the relation (I + h/ {2dfd~ e'‘ , we obtain the approximation 
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2d<{2d{l-f)+{} ~ l + \n(2d{l-f)+l/2. or 

^ < 1 _ 1 +ln(2Q'(l -f)+l/2d. 

2d 


(5) 


Actually, this is a good approximation of 5 which can also be used to show that 
E(«) /, when d and obviously remains d = 0{r). Fortunately, when is in a 

practical range (for example, when the nodes contain from 10 to 400 keys), the expected 
storage utilization is much larger than /, as Table 1 shows. However, it can be noted that 
in this model, storage utilization decreases as the dimension of the nodes increases. These 
results agree with intuitive idea that, if we have r keys {r — > oo), which have to be randomly 
distributed in n cells containing at least 2d/keys and at most 2d keys, the more d increases, 
the more empty space is left in the cells. 

(v) the f= 0.5 case : 

When /=0.5, that is, when the B-tree analogues have their standard formulation, the 
formulae involved simplify and therefore allow us to be more precise. In particular, the 
relation P{t) = 0 can be written in the form 

(1 - = 1 - r or r (1 - r'^) = 1 - 1 ‘‘ 

which simplifies to r (1 + 1 ‘‘) = 1. 

It is worth noting that in the former form, we add the root / = 1 to P{t\ while in the latter, 
we eliminate the di\\ roots of unity. 

Theefore, 1 has the same roots as P{t), except for the di\\ roots of unity. The latter, 

however, are not significant in the asymptotic evaluation of S and W, The real root evaluation 
of i- (1/2 < .s < 1) satisfies the identity 

+ (6) 

and this can be used to simplify all our expressions. We can also easily prove a stronger 
limitation for .s, as follows : 

Using the equation (4) (5) and (6), we get the result 
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Injd^l) In(ln(^+1)) 

(^+1) (^+1) ^ (rf+1) 


(7) 


In fact, if we denote the two values by x and y, i.e. x, s < y, we can show that 

1 - x ~ x^^^ > 0 and 1 ~ y - < 0 which means that x and y separate for root i*. In fact, 

by writing m for ^+1, we find 


In 

_ 

^ Inm^ 

m 

'"in iriy 

- exp 

m In f 1 - 


m 

\ J 


m 

\ J 


m J 

1 1 ' 

[ ~ JJJ 


^In ni^ 


m 
V J 


exp 


( - In m) 


^ (In mf' f In 

2m 

J 


m 


3m" 


+ ... 


^In r 1 ^ 






'^ (In 

2m^ 


+ 0 


(In my 
,2 


m 


since exp (- In /n) = — etc. 

m 


which is positive for m > 3, or d > 2. In a similar way, we can show that 
1 _y_/+i<o. 

By applying the rule = 1 - j, we find the values : 


Q{s) = 


s 

1-s) ’ 


P{s) = 


(2s -1) (d-ds + l) 
s^(l-s} 


and therefore, 


E(u)~ = 


= i2s-l)(d-ds+l) 
2dsil-s) 


1 - 


2s 


1 + 


1 


d(l-s) 
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By means of this formula, we deduce the asymptotic behaviour 
E{u)~ 1/2+ l/21n(cl+ 1), which explains the slow convergence towards 1/2 of the 
expected storage utilization. 

We conclude by pointing out that we only examined the root 5 of minimal modules. As 
previously mentioned, we also have ld-2 complex conjugate roots, some of which have 
modules less than 1. By applying Euclid’s algorithm to P(t) and P'{t), we can show that 
all roots have multiplicity one. These roots produce oscillations in the behaviour of E{u), 
which tends to disappear as r increased. For example, when 10, complex roots having 
minimal modules are : 

s, s 0.831630-0.4155521 

^2 =0.831630-0.4155521 

l^ii = 1 ^ 2 ! H 0.929673 

and the minimal modules real root is 

^3 = 0.8444398 

In this case, 5 and W are more accurately described by the formulae : 


5 ~ 


W ~ 

In which the 
respectively. 


E Ai 
c’ ’ 

/ 


s 

1=1 


Bi{r-2df+ 1) C 


constants correspond to A, By C evaluated in ^ , 


/=1,2,3 
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Abstract 

The motion and stability of a system of two satellites, connected by a light, flexible and inextensible 

string, have been studied under the central gravitational field of the Earth. The non-linear effects due 

to earth oblateness have been analysed using method due to Bogoliubov and Mitropolsky. 

(Keywords : inter connected satellites/earth oblateness/non-linear oscillation) 

Introduction 

The present paper deals with the study of non-linear effects of the Earth’s oblateness, 
on the motion of the system of two satellites connected by a light, flexible and in-extensible 
string in the central gravitational field of force, for the elliptical orbit, relative to the centre 
of mass. 

The relative motion of two satellites connected by a light, flexible and inextensible 
string in the central gravitational field of Earth was studied by Beletsky^’^ and Beletsky 
and Novikova^. The study was based on the assumptions that the orbit of the centre of 
mass of the system is a circular and the system itself moves in the plane of the orbit of the 
centre of mass, Singh^’^ generalized the work by studying it in the two dimensional as well 
as three dimensional cases. 

Manazirudin and Singh^’^, studied the non-linear effects of small external forces on 
the motion and stability of the system of two satellites, connected by a light, flexible and 
in-extensible string. Narayan and Singh^ studied the non-linear oscillation of two 
inter-connected satellites system under the solar radiation pressure. Finally Singh et 
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studied non-hnear parametric resonance oscillation of the system of two inter-connected 
satellites orbiting around an oblate earth. 

We have in this paper, studied the effects of Earth’s oblateness on the non-linear 
oscilition of two artificial satellites connected by a light, flexible and in-extensible string 
relative to its centre of mass moving in elliptical orbit under the central gravitational field 
of oblate Earth. 

A system in space cannot escape perturbing forces and therefore even an infinitesimal 
perturbing force can change the circular orbit into an elliptical one. The complexity of the 
elliptical motion of the system of two inter-connected satellites in space, forced us to restrict 
our analysis in two dimensional case of the motion of the system. 

Equation of Motion 

The equation of two dimensional motion of one of the satellites, when the centre of 
the mass moves along Keplarian elliptical orbit in Nechvill’s coordinate system^^ has been 
obtained in the form Singh 


r' - ir - 3Xp = 


r + 2r = r + 




4AX 


Al 

p 


( 1 ) 


v\ here A'-axis is the direction of the position vector joining the centre of mass and the 
attracting center and T-axis along the direction normal to the position vector in the orbital 
plane of the centre of mass and in the direction of the motion of the satellite. 

In the case the condition for constrained motion is given by the inequality. 

+ F- < A 

P 

1 

1 + ^ cos V 



where 


( 2 ) 
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where p is the focal parameter, e and v are respectively the eccentricity and the true anomaly 
of the orbit of the centre of mass of the system. Prime (/) denotes the differentiation with 

respect to the true anomaly v 


X = p3 A ^ Z! p4 


im^ + m2) 


(3) 


where X is the undetermined Langrange’s multiplier, A is the Earth oblateness parameter 
and jx is the product of mass of attracting centre (Earth) and the gravitional constants. 
When the motion of the satellite m, is determined with the help of eqn. (1), the motion of 
the satellite m 2 is easily determined with the help of identity. 


mi Pi + m2p2 = 0 (4) 

where ^ and ^ are the radius vectors of satellites of masses mi and m 2 respectively with 
respect to the centre of mass of the system. 

As free motion is bound to be converted into constrained motion with lapse of time, 
we shall analyse the constrained motion. In the case of constrained motion the equality 
sign holds in eqn. (2) i.e., the system will be moving along the circle of variable radius 
given by 


= \ (5) 

P 

Now, the equation for motion (1) can be transformed to the polar form, which is best 
suited for our further investigation of the problem given by : 

X = (1 + ^cosv) cost]/ 

7 = (1 + ^cosv) sinxi/ (6) 

where \)/ is the angular deviation of the line joining the satellites with stable position of 
equilibrium of the system. 

Now, differentiating eqn. (6) with respect to true anomaly, putting the derivatives obtained 
above in eqn. (1), multiplying the first equation obtained above by sin \j/ and the second 
equation by cos \(/, subtracting, we get : 
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( 1 + e cosv) \|/" ~ le \|/' sinv + 3 sin\|/ cos\i/ = 2e sinv + 5A (1 + e cosv)^ sin\|/ . cos\|/ (7) 

The eqn. (7) is the equation of motion of a dumbell satellite in the central gravitational 
field of the oblate earth. The eqn. (7) will be the starting point for our subsequent discussion 
of the problem of oscillation. Again, differentiating eqn. (6) with respect to true anomaly 
V, putting the derivatives in eqn. (1) multiplying the first equation obtained above by 
cos\|i and the second equation by sin\|/ and adding we get : 

(i + ecosv)(\j/'+ if-(i -3cosV)--^(^ +£cosv)^(4cosV-sinV) = ^a(^ -f ^?cosv) (8) 

The eqn. (8) determines the undetermined Lagrange’s multiplier The mechanical 
implication of the is that, the motion will be constrained so long as X(r) > 0 i.e. 
X^ (f) < 0. Boundaries between the regions of constrained motion and the regions of free 
motion are given by : 


(w' + 1)3 _ ^ ~ - A {I + ecosv) (4cos V " sin V) = 0 (9) 

(1 + ecosv) 

Equation Characterizing the Non-linear Oscillation of the System 

The non-linear oscillation of the system in the central grvitational field of oblate Earth 
is determined by the equation Singh 


( 1 + ^cosv)\|/'' - leij/' sinv + 3sint|/ cost]/ = 5A (i + ^cosv)^sin\)/ cost]/ + 2e sinv 


( 10 ) 


where 


A = 3 




(Xp 


^8e 


r: 


= oblateness of the Earth = (R^ - R^) / R ^ ; 

Q = angular velocity of the rotation of the Earth; 

= force of the gravity at equator of the Earth; 

Rp = polar radius of the Earth; 

R^ = equatorial radius of the Earth. 

It has been established that there exists a stable position of equilibrium given by : 
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(p = 0 


= 0 


(11) 


The position was found to be Liapunov stable'^. 

In the present paper, assuming e to be a very small parameter instead of zero; Now 
we shall study the oscillation of the system in above mentioned stable position equilibrium. 

¥ = ¥o + B 


¥' 


sinti; = q - ^ + 


( 12 ) 


cost)/ = 1 - 


n! 

21 


Retaining the terms up to third order in the expansion of sin ij/ and cos \|/ in the 
ascending power of \|/, the eqn. (10) takes the form : 

(1 + ^cosv)ri'' - 2^Ti'sinv - 2^sinv - 3ri - 2v^ = 


5Aq - ^ Aq^ + lOAeqcosv - y Aecosvq^ + 5Aqe^ cos^v - y Ae\^ cos^v 


(13) 


Since, the non-linearity has been assumed to be sufficiently weak, 2r\ is supposed to 
be of the order of e. 

Now 2q^ ~ 0(e) and 10/3 Aq^ ~ 0 (e) . 

Setting 2q^ = ek 2 'r\^ ’ ‘"Produce the substitution. 


q = 


(1 •+■ e cosv) 


= (zP) 


( 14 ) 


Differentiating eqn. (14) with respect to v, we get : 
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z" = (l + ecosv)ri" - 2^ri'sinv - ericosv 

z” + eTicosv = (1 + ecosv)ri" - 2eri'sinv (15) 

With the help of eqn. (14) and (15) the equation of motion given by (13) takes the 
form 


+ n~z = e 


-j 20 1 ^ 3 

2z cos V + 5Az cos V + 2 sin v -k^z + — Az cos v - 6z cos v -f ^ 2 ^ 


+ e"' 


- 2z cos^v + 3kiZ^ cos v -f ISAz^cos^v , 


+ e' 




( 16 ) 


Here «“ = 3 - 5A ; where As = (2/e) and = 10/3^ 

The eqn. (16) is the variational equation describing non«linear oscillation of the system 
about the position of equilibrium. We shall use this equation in our further discussions. 

We shall construct an approximate solution of eqn. (16) in the non-resonance case 
exploiting the method^^ of Bogoliubov Krilov and Mitropolsky. After some lengthy 
calculations, we obtain the solution in the first approximation as 

z = u cos e ( 17 ) 

where 0 = nv + 0=5* where, the phase 0* is a constant quantity. The amplitude a and the 
phase 0 ai‘e given by the equations 

^ = 0/.e.(3. = Constant (18) 


d6 3e 
— = n + — 
dv 8n 


(ki 


- k^) ^ 


and in the second approximation the solution is 




4(2«+l) 


4(2n-l) 
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e 


(10A-9)a^ 
12(4/1 -1- l)(2/i+ 1) 


cos(v + 30) + 


( lOA - 9)a^ cos(v - 39)'^ 


(3 - 5A)a^ cos 3el 

12(4/1- 1)(2/!-1) j 

+ 

1 

to 

1 


(19) 


where the amplitude and the phase are determined by the equations : daldv = 0\ a = ag 
(constant) in the second approximation. 


dv 4n 


(4a - IQAg - lOAa^ - 9a^) (12 - 45a^ - IQAa- - 12QA 
96/z (4/1" — 1) 


a^(10QA^-81) 
144(16/1^- l)(4/i^- l)/i 


3e (3 - 5A) 
8 //^ 


( 20 ) 


From the above solution we come to the conclusion that the amplitude a of the system 
remained constant up to the second approximation. The phase of the oscillation of the 
system in the case of non-linear non-resonance oscillation varies with respect to time; 
However, the variation of the phase of the system is of the order of the square of the 
eccentricity, if the orbit of the centre of mass of the system which is a small quantity. 

We arrive at the conclusion that the system experiences main resonance at n = ± 1 the 
parametric resonance at n = ± 1/2 up to the second approximation. 


Non-Linear Oscillation of the System about the Position of Stable Equilibrium 
for Small Eccentricity in the neighborhood of Main Resonance. 


In the previous section, we have observed that the system experiences main resonance 
at n = 1. Now we proceed to construct the asymptotic solution of (16) for the case of main 
resonance n = 1 , in the form : 

c = acos(v + 0) (21) 

where the amplitude a and the phase 0 are given by the system of differential equations : 


da 

dv 


2e 

(«+i) 


cos 0 


dv 


= (n - 1) 4- 


2e 

a{n + 1) 


sin 0“ 


(3 - 5A)a^ 
An 


( 22 ) 
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We conclude that the oscillation of the system in the first approximation is depending 
upon the form of the orbit of the centre of mass. 

The equation can be integrated in a close form as it can be reduced to the Hamilton 
canonical form; 


where 


da 

dv 


a 


dH'] . dQ 

36 = 

V J 


I 

a 


(M] 

da 

V J 


H = -2a 


(n+l) 


• A (n-l) 2 3-5/1 4 

sin 0 - — ^ -2 + -^- 77 — - ^ 

2 16n 


(23) 

(24) 


Obviously, the system of the eqn. (23) has a first integral of the form H = Co (25) 

which reduces the problem to quadrature. Here Co' is the constant of integration. However, 
it is preferable to analyse the integral curves in the phase plane (a,Q). In order to plot the 
integral curves in the phase plane (a,0) let us put the eqn. (24) in the form : 


in-- 1)/ + 4aesin0--^^^^(3-5/l)a^ + Co = O 

o/Z 


(26) 


where Co = 2 (n + l) Co' 

The integral curve (26) has been plotted in Fig. 1 for « = 0.8, e = 0.01 and /\ = 0.18 
(i.e. Earth oblateness parameter). The integral curves drawn in the phase plane (a, 0) clearly 
indicate that there exists only one stationary regime of the amplitude i.e. a = 0.046 and 0 
= (81)° and it is stable as the integral curves are closed curves. For any other initial 
conditions, we shall obtain periodic change in the amplitude a which will be bounded. But 
the maximum values of a in this case will always be greater than its values at the stationary 
regime. We conclude therefore, that for a gravity gradient stabilisation of such a space 
system m elliptical orbit, we need to bring the amplitude of oscillation near the stationry 
regime which gives, smallest deflection of the system from the relative equilibrium position 
in comparison, to any other regime of oscillation. 

In the Fig. 2 the integral curves, (26) for «= 1 .2, e = 0.01 and A = 0. 1 8 (Earth oblateness 
parameters) has been plotted and it is found that the stationary regime of the amplitude a 
exists. In this case also there exists only one stationary regime which gives, smallest 
e ection of the system from the relative equilibrium position in comparison, to any other 
regime of oscillation. We also conclude that the period passes through the position of main 
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resonance stationary amplitude declines steadily, which is obvious from Fig. ( 1 & 2) plotted 
for n = 0.8 and n- 1.2 respectively. 

In order to see, whether the oscillation of the system in this stationary regime always 
takes place with tight string, we shall consider the eqn. (9). This condition divides the phase 
plane into two regions which are regions for free and constrained motions. The boundary 
of these regions is given by : 

(1 + ecosv) {r\' + 1)^ ~ Scos^T] + 1 - A (1 + ^cosv)^ (5 cos^r) - 1) = 0 (27) 

cos'll = cos^ (acQsk) = cos^ [acos (v + 6)] 

Jq (a) - 27, (a) cos(20 + 2v) + 2J^ (a) cos (40 + 4v) 

2 2 

l-y-YCOs(26 + 2v) 

sin Tj = 27^ (a) cos k - 2J^ (a) cos 3k + 27^ (a) cos 5k (28) 

= cos (v + 0) 


where 




(- D" 


/ V« + Y 

' a ' 


m=0 


m 


:Vy+W+ 1 


(29) 


Now substituting the value of T|, r\' cos'r] in the above eqn. (25) neglecting the terms 
containing e" and higher order of e and averaging with respect to v in the interval (0, 271 ) 
keeping 9 and A to be constant we get, 


aV (3-5A) 
2 4 


ane 


p- COS0 - 1 + ^ (2^ - 4/4 + 

(n + 1) 2 2 


- ane sin i 


= 0 


(30) 


In the process of tracing out the average eqn. (30) for n = 0.8, e = 0.01 and A = 0.18 
find that the boundary line does not intersect any of the integral curves but passes just above 
them. Also, from the figure in which the averaged equation (30) for n = 1.2, e = 0.01 and 
A = 0.18 has been traced we find that the line does not intersect the integral curves. 

Therefore, the motion of the system in the whole phase plane (a, 0) will be 
non-evolutional. 
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Discussion and Conclusion 

We thus, conclude that the amplitude of the oscillation of the system in the case of 
non-linear non-resonance oscillation remained constant up to the second approximation. 
The phase of the oscillation of the system varies with respect to time. However, the variation 
of the phase of the system is of the order of the square of the eccentricity, if the orbit of 
the centre of mass of the system which is a small quantity. We observe that the system 
experiences main resonance at /2 = ± 1 and parametric resonance at n = ± 1/2 up to the 
second approximation. 

We further conclude that if a system of two cable connected satellites oscillating at 
and near the main resonance n = 1 under the gravitational field of oblate Earth, there exists 
only one stationary stable regime of amplitude for both /2 < 1 and/ 2 > 1. We have established 
that at the resonance oscillation at and near 22 = 1, the motion remains non-evolutional i.e. 
the system moves like dumbell satellite. However, the stationary amplitude declines steadily 
as it passes through two values of n given by w < 1 and 22 > 1. Thus, we arrive at the 
conclusion that oblateness of the Earth, however small, may cause periodic changes in the 
amplitude of the oscillating system in case of slight deflection from stationary regime of 
oscillation. 
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Abstract 


A unique common fixed point theorem for a sequence of selfmaps of a complete metric space by 
altering distances between the points under a certain continuous control function is established. 

(Keywords : complete metric space/Hilbert space/selfmaps/common fixed points) 


Introduction 

Existence and uniqueness of fixed points for selfmaps of a metric space by altering 
distances between the points with the use of a certain control function was established by 
Khan et al.^ Some contributions in this direction are due to Pathak and Sharma“, Sastry 
and Babu^’"^. 

Throughout this paper, (X, d) denotes a complete metric space, H a Hilbert space, 
the set of all non-negative real numbers and N the set of all positive integers. <& denotes 
the set of all continuous selfmaps ([) of satisfying (i) ^ is increasing (ii) <j) (0 = 0 if and 
only if r = 0. 

Existence and uniqueness of fixed points for a sequence of selfmaps of a Hilbert space 
H is established by Pandhare and Waghmode^. 

f I®® 

Theorem 1 : Pandhare and Waghmode.^ : Let C be a closed subset of H and 7; 

I >n—\ 

a sequence of selfmaps of C satisfying 
(A) : For any two maps , Tj 

\\T.x - Tjyf < alk-yip + b(\\x-T^xf + ||y-r^|P) 
for all distinct x, y in C where a^O, 0 < b < 1 and a-^2b < 1 . 
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(B) : There is a point Xq e C such that any two consecutive members of the sequence 
{a-„} defined by , n > 1 are distinct. 

I I cso 

has a unique common fixed point in C. 

1^=1 

Babu^ obtained a partial generalization of theorem 1 to complete metric spaces by 
assuming more stringent hypothesis that all the members of the sequence in (B) are 
distinct. In this paper, we obtain a generalization to theorem 1 in complete metric spaces 
under hypothesis (Bl) (stated in hypothesis of theorem 2) by altering distinces between 
the points, with the use of a certain continuous control function ^ in <I), which is our main 
theorem. 

Examples are also provided to justify the study of the mappings under condition (A) 
(in the Hilbert space setting) and the generalized version under condition (Al) in metric 
space setting. In fact, Example 3 shows the applicability of theorem 2 in /|-norm setting. 

Main Theorem 

Theorem 2 : Let \tT be a sequence of self maps on (A, d). Assume 

(Al) : There exists a <}) in O such that 

(t)(d (T; x, Tj y)) < a (t) {d (jt, y)) + b {^{d {x, 7,. x)) + (?{d{y, Tj y))) 
for all ij in N and for all distinct x, y in X where a > 0, 0 < Z? < i 
with (3 + 2Z?< 1. 

(Bl) : There is a point Xq in X such that any two consecutive members of the sequence 
defined by , n > 1 are distinct. 

Then |T„| has a unique common fixed point in X. Infact, |x,^} is Cauchy and the limit 
of lx,] is the unique common fixed point of jr,,} 

Proof : Write a„ = d (x„,x„^i) and P„ = (f)(a„). From (Al) and (Bl), we have 
^i = ^{d{x„X2)) = i?idiTiXQ,T^0 


< a (j) (of (xq , xj)) + b(^id (xq, x,)) + (f» (J (x,, Xj))) 
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so that (3[ < k (3q where k = < I, since a + 2b<l . 

By induction, follows that ^ ^ P„_i for all n > 1 (1) 

so that i 0 as 7 z °o and a„ < a„_| for n = 1, 2, ... Therefore {a„} is a decreasing 
sequence of non-negative reals. Thus {a„} i a (say). Then 

(ct) so that (t)(a) = 0 and hence a = 0. Therefoe {a„} 4- 0. (2) 

The crucial part of the theorem is to show that {x„} is Cauchy in X. Otherwise, there 
is an e >0 and sequences [m (k)} and { « (^) } such that m{k) < n(k), d ^ 

Assume that =-r, „(/.)_] for infinitely many k. Then, for such k, 

we have 

€ < < d + d x, „(;,)) 

= ^ i^k)’ ^n(q-i) + d 0 as * ^ oo , a contradiction. 

Hence, for large fc, A , consequently 

= <t> (^n('t) ^n(k)-V '^m{,k) ^m(k)-\)) 

< a<i){d{x„^k)-x )) + ^ (^ i^n{k)-V ^n{k^) + 

< <2 (j) (e -H ^,„(i)-i)) + b i.d + ^{d x^(i)))) 

a (j) (e ) as A: -> o°, by (2). 

Hence (|)(g ) < a ([)(£ ) < (J)(e ), a contradiction. 

This shows that {x^} is a Cauchy sequence in X. As X is complete, limit of {x„ ) exists. 
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There is a sequence {w(^)} such that y • Otherwise, y = for large n, which 

is not the case, since consecutive terms are different. With this subsequence we 

have, for any positive integer m, 

= ( 1 ) 

<a(^id(y, + bi^{d {T^, y)) + (j) (d Jc„(*)-i))) 

Taking limits as k^°°, we have (|) (d (T^, y))<b^ {d {T^, y)). 

Since 0 < Z? < 1 , it follows that <[) {d {T^, y)) = 0 so that d {T„y, y) = 0. 

This shows that y is a fixed point for . Thus >> is a common fixed point for the sequence 

{7’„ir=i • 

Uniqueness of the fixed point follows trivially from (Al). 

Remark : When X is a Hilbert space, we get theorem 1 due to Pandhare and Waghmode^ 
from the above theorem, by taking ^{t) = for r > 0. 

Generalization of Theorem 2 

Theorem 2 can be further generalized (the proof being similar) as follows. 

Theorem 3 : Let be a sequence of selfmaps on (X, d) and assume (Bl) holds. 

Further assume that 

(A2) : There exists (j) in O such that 

<^{d {Tpc, Tjy)) < k max {^(d (x, y)), ^ ((j) (d (x, Tfc)) + <S^(d(y, Tjy ))) } 

for some 0 < Z: < 1, and for all ij in N and for all distinct x, y in X. 

Then the sequence {7^}“=! has a unique common fixed point in X. In fact, {?„} {x„} 
is Cauchy and the limit of {TJ {x„) is the unique common fixed point of 

Examples 

The following is an example to show the applicability of theorem 2 with (l)(r) = t. 
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Example I ; Let X = [0, 0.2] with the usual metric. Define T„ : X by 

for « = 1, 2, ... Define (j) (r) = r, r>0 that <j)G <5. Then, {TJ” , satisfies the condition 
(Al) with r/ = 0.5 and 6 = 0.125. Observe that, for any non-zero Xq in X, the sequence 

{•'^,,1,7=1 defined by x„ = n > 1 has all its elements distinct so that (Bl) also holds; 

thus hypothesis of theorem 2 is satisfied and 0 is the unique common fixed point of 

(7'J7=i- 

Example 2, below shows the applicability of theorem 2 with ^(t) = p-. This example 
may be taken as an application of theorem I in Hilbert space setting. 

Example 2 : Let X = [0, 0. 1] with usual metric. Define : X X by T,p: = for 
11 = 1,2.... Define (]) (/) = P, r > 0 so that (|) G <I). Let x, y g X, x^y . Then 

0 (d (T,^, T,„y)) = (x^" - /"f < (0.04) ex'' - ff < (0.04) (x^" + /") (3) 

and (0.04) (x^'' +y^"') + 2 (0.05) (x^"^‘ + y^^') 

= (0.04 + (0.1)x)x2"-t-(0.04 + (0.1)y)/'" 

< (0.05) (x^" + y^"') 

< (0.05) (x^ + y^) 

< (0.01) (x-y)^ -1- (0.05) (x^ y^) + (0.05) (x'*" -t y'*'”) 

for all m > 1 and n > 1, 

so that (0.04)(x^'' + y^'") < (0.01) (x-y)^ + (0.05) [(x-x^^ + (4) 

From (3) and (4), it follows that the inequality (Al) holds with a = 0.01 and b = 0.05. 

Condition (Bl) holds trivially for any 0 t^Xq G X, and 0 is the unique common fixed 
point for 

In the above example, we observe that (Al) does not hold with (t)(r) = t and a = 0.01 
and b = 0.05 (take ;c = 0 and }^ = 0.1). 

We now give an example as an application of theorem 2, with (})(/) = ? in metric space 
setting, in general and /j-norm setting, in particular. 
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Example 3 ; Let X = 




' 1 ' 

r 

X 

1 

18 


with /[-norm, || . ||j being defined by 


Define T„\X-^Xhy {x, y) = (/", /") for n = 1 , 2, . . . Then { T„ } satisfies (A 1 ) with 
a = 0.5 and b - 0.1. Condition (Bl) holds trivially for any nonzero (.rQ, y^) in X and (0, 0) 
is the unique common fixed point for 


Open Problem 

Is the conclusion of theorem 2 valid if (Al) is replaced by (A3) : There exists 0 in 
<l> such that 

0 [d {Tfc, Tjy)) <a^(d (x, y)) + Z? (0 (J (x, Tpc)) + 0 (J (y, T-y))) 

■f c (0 (d (x, Tjy)) + 0 (li (y, Tp:))) 

for all /, j in N and for all distinct x, y in X where ^>0, c>0, 0<Z?<i with 
a + 2b 2c < 1. 
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Abstract 


In the present paper Hall current effects on unsteady free convection and mass transfer flow of a 
VISCOUS, electrically conducting fluid past an infinite, porous, non-conducting, vertical plate under 
the intluence of uniform transverse magnetic field in the presence of constant normal suction velocity 
at the plate is considered. It is assumed that the free stream, plate temperature, species concentration 
osci Hates about a non-zero mean and the concentration level of the diffusing species is small compared 
with other concentrations. Momentum, energy and concentration equations are solved. The effects 
of the various parameters, entering into the problem, on the primary and secondary velocities are 
depicted giaphically while the skin-friction coefficient, Nusselt number and Sherwood number are 
shown in tabular form followed by a discussion. 

(Keywords : Hall current/mass transfer/heat source/Soret number/Sherwood number) 


Introduction 

The studies on MHD free convection analysis appear in literature because of its 
important applications in meteorology, solar physics, cosmic fluid dynamics, motion of 
earths core, planetary magnetospheres, aeronautics, chemical engineering, electronics etc. 
Along with the free convection currents caused by the temperature difference, the flow is 
also affected by the difference in concentration on material constitution. Gebhart and Pera 
made extensive studies of combined heat and mass transfer flow. The effect of suction 
parameter on free convection and mass transfer MHD flow being an effective method of 
controlling the boundary layer has been studied by several authors . In all these studies 
thermal effect has been neglected. Eckert and Drake^ have pointed out that in a convective 
fluid when the flow of mass is caused by a temperature difference one can not neglect the 
thermal diffusion effect (commonly known as Soret effect) due to its practical application 
in engineering sciences. For instance thermal diffusion effect has been utilized for isotope 
separation and in mixtures between gases with very light molecular weight (H2, He) and 
medium molecular weight (N2, air) and it was found to be of a such magnitude that it can 
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not be neglected. In view of the importance of this diffusion thermo effect, Jha and Singh^ 
presented study for a free convection and mass transfer flow past an infinite vertical plate 
moving impulsively in its own plane taking into account the Soret effect. Kafoussias^^ 
extended this work to MHD free convective flow. Recently Sattar and Alam^ ^ have discussed 
Soret effects on MHD free convection and mass transfer flow past an impulsively started 
vertical porous plate in a rotating fluid with out taking into account the Hall current effect. 
However, when the strength of the magnetic field is strong, the Hail effects play a significant 
role in determining the flow feature. More recently Singh and Singh^“ have discussed Hall 
current effects in convective flow of a viscous liquid. Hence our aim is to study Soret 
effects as well as heat source effects on MHD free convection and mass transform flow of 
a viscous fluid past an infinite porous non-conducting vertical plate with constant suction 
normal to the plate taking Hall currents into account. 


Formulation of the Problem 


Let .r-axis be chosen along the plate and y-axis normal to it. A uniform magnetic field 
Bq is taken to be acting along the y-axis and the induced magnetic field is assumed to be 
negligible so that B = (0, Hq, 0). The plate temperature and concentration oscillate about 
a constant mean and the species concentration is assumed at low level. The plate temperature 
and concentration are instantaneously raised or lowered to (> T^) and Q, (> C^), 
and are the temperature and concentration of the uniform flow. The equation of 

conservation of electric charge is = 0 which gives = constant, where 
— > 

J = {J^, Jy, J.). Since the plate is electrically non-conducting, 7^, = 0 and is zero every where 
in the flow. When the magnetic field is large, the generalized Ohms law, in absence of 
electric field, neglecting the ion slip and thermo electric effect (Meyer'^) yields 


J. - 


( 1 ) 

+ 


( 2 ) 


The solution of eqn. (1) and (2) is 


oix.Hq 

= --- I {mu - vv) 


1 + 

1 + 


{u + mw) 


(3) 


(4) 
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where m{= zj is Hall parameter, : the cyclotron frequency, z^ : the electron collision 
time, [i, : the magnetic permeability, a : the electrical conductivity of the fluid. Within the 
frame work of such assumptions, neglecting the Joule heating and viscous dissipation terms 
and assuming that all the fluid properties are constant except that the density variation with 
temperature is considered only in the body force term, the basic equations governing the 
boundary layer flow past the plate under Boussinesq’s approximation in non-dimensional 


(ignoring 

the stars 

introduced ahead) are 




J. 

dll 

du 

1 du d^u M 

(m- 

- U + mw) + r + 6 


4 

dt 

dy ~ 

4 dt 3^2 1 ^2 

(5) 

1 

dw 

dw 

d^w M _ , 

U)- 

- w] 


4 

dt 


2 + , 2 t'” ~ 

ay 1 + m 

(6) 

ji 

dT 


1 d^T „„ 




4 

dt 

dy 

Pr dy^ ^ 



(7) 

1 

dd 

ae 

1 a^e dh 



(8) 

4 

dt 

ay = 

Sc ay' “ ay' 




The non-dimensional boundary conditions are 

u = w = Q, r=9=l + 6 exp(mr) at y = 0 

M = 1 + 6 exp(mr), vv-^O t = Q = 0 as 


(9) 


where : Grashof number, : modified Grashof number, : Prandtl number, : 
Schmidt number, Sq : Soret number, S : heat source parameter, P : volumetric coefficient 

of thermal expansion, P* : volumetric coefficient of thermal expansion with concentration, 
K : thermal conductivity, D : molecular diffusivity, Dj : thermal diffusivity and all the 
other quantities have their usual meanings' 

The non-dimensional quantities introduced in the above equations are defined as 


» U * W * rf — ti* - 

IT ' Ti ^ y A, ’ 4u ’ [A -y 


Uo' Uo 
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t = 


M 


T-T„ 

G-C^ 

a* ^ o® 


T.-TJ 

n ^ r* ’ 

'^W '^oo 

= - 

VqP 

S - ^ P 
^0 

li 






. G„ 


Govl 


D 




Introducing q = u + iw, the equations (5) and (6) transform 


to 


^ M(1 -im) ^ 

4 3t dy 4 dt 3 / (1 + ^2^ G) + G;r G„Q ( 10 ) 


Solution of the Problem 


In order to solve the system of eqn. (7), (8) and (10) under 
(9), we assume 


the boundary conditions 


<?(>'> 0 = (1 -%) + & (1 -q^) exp(/m) 
7’(y, 0 = TgO) + 6 Tj (j7) exp(m/) 

6(y, t) = 6o (y) + 6 6, 0) exp (int) 


( 11 ) 


Substituting (1 1) in (7), (8), (10) and equating harmonic and non-harmnmV 

get six equations. The solutions of these equations under the bonnHarv a - ■ 

after substituting (1 1) in (9) give conditions obtained 

St Dm give % . . /q - ^ 1 . 00 and Si • Usmg these in (1 1), we get 


q = {M^- i'M,.) + G (I^ - /I.) exp {int) 

T = exp (-Coy) + € (X, - /K,) exp (mr) 
e = (1 4- 4 W exp (- 5,y) - exp (- Coy) -f g (iV, - i;v,) exp (/«,) 
where M,= l-D, + a, (G, - GJoS,) exp (- Coy) + Ug exp (- 5,y) 
^3 ^4 (Gr - G^SqS^) exp (- Coy) + b^G^ exp (- Sj) 


( 12 ) 

(13) 


(14) 
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L, = 1 - D4 + G, Cl + (1 + d,SoS,) G^Cj - £3 G,„SoS, exp (- Coy) 

L. = F4 + Gp^ + (1 + cfiTO G^Dj - exp (- Coy) 

Ni- = X 2 + iX2~ X^) CjSqS^ + e4So5j. (y2 ” Ej) 

N, = ^2 + (^2 ~ E]) e3So5^ — ^4SoSc{X2~ X^) 

{X„ , E„) = (cos y , sin |3^) exp (- a,^) n = L 2 

(C„ , D„) = (A„ cos P„ y + S„ sin p^) exp (- a,j;) n = 1, 2, 3, 4 

F„ = {B„ cos p„y - A„ sin p„ y) exp (- a„y) n = 3, 4 

a„ = al- h ” = 

^. = 2a„p„ - P„ + £2 
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aj 63 

1^3 j ^4 *” 2 2 y ^4 ^ ^ < 

^3 + I>3 413 + Z>3 




_ {l+d,S^,)E, 

*“ 9 ■“ 2 2 

£l + i)3 


(H-di 5 (^,)I >3 


_ y. ■ -3 P _ c2 _ C t, 


^2 ~ ^3 ~ 4 > £3 - 2 ,c 2 ' ^4 ~ 2 , p2 ’ ^1 ~ Pi 

^ ^3 -^2 ^3 "^2 


^2 — 2CX|(3j 5^Pl + 


(a^ - 13^)6, + 2^2 a, Pi 


g] + ^2 


2«i(XiPi - (cti - Pi) ^2 


+ fi 2 
^1 ^ ^2 


A3 - 1 - 04 + G„ {SfjS^a^ - flg) 


B3 = GA - G„ (SoS^5 -b,), A4 = 1 - AiG, - A2G„ - G„So 5 , (^1^2 - £3) 

B4 = Bi G^ + Bj <^m + (^1^2 + ^4) 

Shear Stress, Nusselt Number and Sherwood Number 

The non-dimensional expressions for shear stress x, rate of heat transfer in terms of 
Nusselt number Nu and rate of mass transfer in terms of Sherwood number Sh are 


Nu = - 
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Table I - Transient primary (Xp) and secondary (x^) shear stress. 



(€ =0.0i,n 

o 

II 

= 0.71 and ni 

= 71/2) 





m 

M 

S 

Gr 

G„ 

So 

Sc 

Xp 

Xs 

0.5 

5.0 

1.0 

5.0 

4.0 

1.0 

0.60 

1.4695 

0.6937 

1.0 

5.0 

1.0 

5.0 

4.0 

1.0 

0.60 

1.3932 

0.7528 

0.5 

10.0 

1.0 

5.0 

4.0 

1.0 

0.60 

1.3887 

0.5167 

0.5 

5.0 

2.0 

5.0 

4.0 

1.0 

0.60 

1.5482 

0.6521 

0.5 

5.0 

-2.0 

5.0 

4.0 

1.0 

0.60 

1.3817 

0.7289 

0.5 

5.0 

1.0 

iO.O 

4.0 

1.0 

0.60 

1.3548 

0.5098 

0.5 

5.0 

1.0 

-5.0 

4.0 

1.0 

0.60 

1.5897 

0.7139 

0.5 

5.0 

1.0 

5.0 

8.0 

1.0 

0.60 

1.3482 

0.4973 

0.5 

5.0 

1.0 

5.0 

4.0 

2.0 

0.60 

1.3466 

0.4836 

0.5 

5.0 

1.0 

5.0 

4.0 

1.0 

0.75 

1.6183 

0.6892 


Table 2 - Primary {Nup) and secondary {Nus) Nusselt number. 
(€ =0.01 and nr = 71/2) 


S 

Pr 

n 

Nup 

Nus 

0.10 

0.71 

5.0 

0.9236 

1.0422 

0.15 

0.71 

5.0 

0.7815 

1.0296 

-0.10 

0.71 

5.0 

0.4921 

1.1018 

1.00 

7.00 

5.0 

27.9842 

6.2801 

0.10 

7.00 

5.0 

53.1879 

7.1088 

0.10 

0.71 

10.0 

0.9208 

1.3118 

0.10 

0.71 

15.0 

0.9187 

1.5209 
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Table 3 - Transient primary (Shp) and secondary (Shs) Sherwood number. 
(6 =0.01, Pr = 0.71,n = 5.0 and nr = Ji/2) 


5 o 

Sc 

S 

o 

X 

10 '^ X Shs 

0.0 

0.22 

0.10 

0.3358 

0.5196 

1.0 

0.22 

0.10 

0.2895 

0.4206 

2.0 

0.22 

0.10 

0.2433 

0.3216 

0.0 

0.60 

0.10 

0.5548 

0.9764 

1.0 

0.60 

0.10 

0.5384 

0.9723 

2.0 

0.60 

0.10 

0.5331 

0.9681 

0.0 

0.75 

0.10 

0.6199 

1.1313 

1.0 

0.75 

0.10 

0.5621 

1.2055 

2.0 

0.75 

0.10 

0.5042 

1.2796 

i.O 

0.22 

0.15 

0.2855 

0.4001 

1.0 

0.22 

- 0.10 

0.3001 

0.3843 

1.0 

0.22 

- 0.15 

0.3053 

0.3819 


Discussion 

To understand the physical depth of the problem, some numerical calculations are 
carried out for the non-dimensional primary velocity field (u) and secondary velocity field 
(w). To be realistic, the values of the Prandtl number (P,.) is chosen to be 0.7 1 and 7.00 
which corresponds to water vapour and water respectively. The value of Schimdt number 
(S^) is chosen to represent hydrogen {S^ = 0.22), water vapour (S^ = 0.60) and oxygen 
(S^ = 0.75) respectively. The values of Grashof number (Gr), modified Grashof number 
(Gm), Soret number (Sq), magnetic parameter (M), Hall parameter (m), heat source parameter 
(S) and frequency parameter (n) are chosen arbitrarily. 

The variations of primary velocity (k) and secondary velocity (w) are shown in Fig. 
(1&2) respectively. A deep study of these figures reveals that (i) an increase in M leads to 
a decrease in primary velocity but an increase in secondary velocity, (ii) an increase in 
Soret number Sq increases the primary velocity but decreases the secondary velocity, (iii) 
an increase in heat source parameter (absorption type) increases primary velocity but 
"decreases the secondary velocity and vice-versa for generation type heat source parameter. 
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(P^= 0.71, (3^= 5.0, = 4.0, n = 10.0, and 8 = 0.01) 


Curve S M So 



m S(. 
0.5 0.22 
0.5 0.22 
0.5 0.22 
0.5 0.22 
0.5 0.22 
1.0 0.22 
0.5 0.60 


- 0.71, G^= 5.0, = 4.0, n = 10.0, and 8 = 0.01) 


-0.05 


- 0.10 



-0.15h 


-0.20 h 


-0.25 


Curve S 

M 

s. 

m 

Sc 

1 

1.0 

5.0 

1.0 

0.5 

0.22 

II 

2.0 

5.0 

1.0 

0.5 

0.22 

III 

-1.0 

5.0 

1.0 

0.5 

0.22 

IV 

1.0 

10.0 

1.0 

0.5 

0.22 

V 

1.0 

5.0 

2.0 

0.5 

0.22 

VI 

1.0 

5.0 

1.0 

1.0 

0.22 

VII 

1.0 

5.0 

1.0 

0.5 

0.60 


4.0 


Fig. 2 - Secondary velocity profiles. 
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(iv) an increase in m increases primary as well as secondary velocity, (v) an increase in 
decrease primary velocity but increases the secondary velocity. 

The variations in non-dimensional transient primary and secondary shear-stress 
(denoted by and x^), rate of primary and secondary heat transfer in terms of Nusselt 
number (denoted by Nu^ , transient primary and secondary mass transfer rate (denoted 
by Shp and Sh^ are shown in Tables 1, 2 and 3 respectively. The conclusion regarding the 
behaviour of the parameters on shear-stress, Nusselt number and the Sherwood number 
are self evident from the tables. Therefore any further discussion about them seems to be 
redundant. 
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Abstract 

The aim of this paper is to investigate the flow of visco-elastic Oidroyd fluid with transient pressure 
gradient through a rectangular channel. Firstly, the first order Oidroyd fluid has been considered. 
Secondly, the authors propose a second order Oidroyd fluid and investigate the same problem in this 
new visco-elasiic model. The numerical computation of the velocity profiles have been presented for 
both the first and second order fluids in the form of tables and graphs. The result for visco-elastic 
Maxwell fluid and purely viscous fluid have been determined from the first order Oidroyd fluid simply 
by passing to the limit as p, 0 (Maxwell fluid) and A, 0, p, 0 (purely viscous fluid). Finally 
a new generalized model of visco-elastic fluid has been proposed. 

(Keywords : visco-elastic Oidroyd fluid/unsteady flow/transient pressure gradient/rectangular 
channel/generalized visco-elastic fluid) 

Introduction 

In the classical theory of viscous flow tremendous development have been observed 
as evidenced by the informative monographs of Lamb ^ and various other authors^ In 
recent years the development of new areas of fluid mechanics are very remarkable. For 
this one may refer to the review literature^ in connection with non-Newtonian fluids, 
polymeric liquids and visco-elastic liquids of different types. Moreover some interesting 
problems in this area have been investigated by Sengupta and his research collaborators 

Basic Theory and Equations of Motions for First Order Oidroyd Fluid 

For slow motion, the rheological equations for Oidroyd visco-elastic fluid are 


= - p'h + 
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1 + Xj 

V 





2|i 


l+|ii 



) 







where is the stress tensor, is the deviatoric stress tensor, eij the rate of strain tensor, 
p' the pressure, X/ the stress relaxation time parameter, jLi/ the strain rate retardation time 
parameter, 5^- the metric tensor in cartesian co-ordinates and |i, the co-efficient of viscosity 
and V- the velocity components. 

Let us now consider the walls of the rectangular channel to be the planes x = ± a 
and y = ± b where z'-axis is taken towards the direction of motion. 0, 0, w' (x\ y/) are 
respectively the velocity components along x\ y\ z directions where w' (/, y/) is axial 
velocity of the fluid. A transient pressure gradient - Pe~^^ varying with time is applied 
to the fluid. 


Following the stress-strain relations, the equation for unsteady motion is given by 

,3^ 


1 + 


3 /' 


8w' 

dt' 


I 

P 


l+Ai 






V^w' 


( 1 ) 


Introducing the non-dimensional quantities 


vv = vv 


/ a 


P^P ~ 
pv^ 




“ , 0) = CO — , (x, y, z) = - (x , y , z ), 





, V 

a 


eqn. (1) becomes 


1 -f Xi 


dt 


dw 
dt '' 


1 + X| 


dt 


dz 




w 


(2) 
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We have considered those type of situations of the flow which is transient in nature 
with respect to time and periodic in nature with respect toy. From the nature of the boundary 
conditions we have chosen the solution of (2) as 


w - W (x) cos mye~^^ 

(3) 

The boundary conditions of the fluid are 


(i) W = 0 when x = ± 1, - ~ < y < ~ 

a a 

(4) 

b 

(ii) w = 0 when y = ±— ,-~l<x<l 

a 

(5) 

Condition (5) will be satisfied if cos ??z - = 0 

a 


or m = {2n -hi) ^ n = 0, 1, 2, 3 

lb 

(6) 


We construct the solution as the sum of all possible solutions for each value of n. 


vv = ^ W(;c) cos mye'^^ 


(7) 


By putting 


dz 


Pe~^^ ( ol ) > 0) in (2) and using (7), we get, 


Z 


Wix) 2 „„ ,1 
— — W(.x)>cosmy + 


oo 

I 


co(l -X, (o) 

(1 - H,(o) 


W{x) cos my + 


P(l-Xt(0) 

(l-JilCO) 


= 0 


and equating the co-efficient of cos my equal to zero, the value of W can be determined 
from 


dx^ 




VV + A„ = 0 


(9) 
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A„ = 


(-1)" 4P(1-A,,(0) 


where 


(2n+l) n(l-iJ,, CO) 


= W 


(0(1 -X,(0) 


(1-Pl®)j 

Using the boundary condition (4), the solution is 


W{x) = 


(- 1)" 4P(1 -X,co) 


(2n+l) Ji(l-p,i(o)/!:^ 
So the velocity of the fluid is given by 
(-1)" 4F(l-XiC0) 


1 


u ^ 

cosh — X 
a 


cosh 


K 


W{x,y,t) = X 
;?=0 


1 ~ 


cosh ^ 
a 


cosh — 
a 


Basic Theory and Equations of Motions for Second Order Oldroyd Fluid 


e-""' cos(2n+l) 


nay 

2b 


(10) 


( 11 ) 


For slow motion, the rheological equations for second order Oldroyd visco-elastic fluid 


are 


T.. = -p' 8 y + T/ 

1 +Vt7 + V-^ 
‘ dl' 2 3^,2 


V = 2ii 


1 + 1^1 37 + P2 


, a- 


,2 ^ 


dt' ■ ^ 3/2 


■ij = 0 ('’•J + 


where Xy is the stress tensor, Xy is the deviatoiic stress tensor, ey the rate of strain tensor, 
p' the pressure, the stress relaxation time parameter, the strain rate retardation time 
parameter, Xj' is the additional material constant, is the additional material constant, 
5y the metric tensor in cartesian co-ordinates and ft, the co-efficient of viscosity and v,- the 
velocity components. 

Let us now consider the walls of the rectangular channel to be the planes x' = ± a 
and y' = ±b where z'-axis is taken towards the direction of motion. 0, 0, w' {/, y', t') are 
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respectively the velocity components along x\ y\ z directions where w' (x', y', t') is axial 
velocity of the fluid. A transient pressure gradient - Pe~^‘ varying with time is applied to 
the fluid. 

Following the stress-strain relations, the equation for unsteady motion is given by 


1 + A. 


at' ^ dt'y dt' 





'^1 

y 



+ V 


l + Fi' 


a , a-^ 


V^w' 


(12) 


Introducing the non-dimensional quantities in eqn. (12) 


2 2 

/Cl / G' /V / G . \ l/ ^/^\ 

w = w - , p = p — r , t = t , 0) = (0 — , {x,y,z) = - {x ,y ,z), 
V G^ V a 


w V , V T 

~2’ Fl ~ Fl 2’ ^2 -^2 4’ I*’2~I*2 4 

a a Cl a 


eqn. (12) becomes 


42 ^ 


1 -f A. 1 -V + X9 




dw 

dt 


l+^l^ + ^2^ 


+ v 


+ ^^2 




VV (13) 


Solution of the Problem 


Let us assume the solution of (13) as 


w = W{x) cos mye~^^ 

Here we have used the same boundary conditions of the fluid given by (4) and (5) and 
used the relation (6). 

We construct the solution as the sum of all possible solutions for each value of n. 
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vv 


= ^ W(x) cos my e~ 

;i=0 


(14) 


Using (14) and putting 


ino - pe'*^ (to > 0) in (13), we get, 


E 

/i=0 


cr W{x) 2 irr/ N V 

^ - m W(x) f cos my + Ji, 

dX' 


ca(l -XiC0 + X2C0^) 


. W{x) cos my + 

„=0 (1 ) 


If we express 


P(1 — A,j(jl) + ^2^^) 
(I - p,a) + lajco^) 

P (I — X|C0 + X2tu^) 

(i — Xj(l) + X20)*^) 


(15) 


as a Fourier cosine series in the interval of 


- “ < V < - and equate the co-efficient of cos my to zero, the value of W can be determined 
a ' a 

^ d"W a 

from — 7 ~ — W + = 0 


dx^ 


where 






lY 4P (1 - XjO) + ^2^^) 
(2n + l) 7i(l -|ij(o + |i ,2 00^) 




(0(1 ~Xia) + X2a)^) 
(I - 1 X 1 ( 0 x 1 X 20 )^) 


We now obtain the solution by using the boundary condition (4) as follows 

W(x) 


(~ 1)^ 4P(i -X j(0 + A,2C0^) 

(2n +1) 7t (1 - 1 x ^(0 -H |i 2 ( 0 ^) 


1 - cosh™ X 
a 


cosh 


K 


So the velocity of the fluid is given by 
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J.O = X 

n=0 


(■„ lY 4P (1 - XiCOH- 
(2n + 1) 71 (i ~ ii|CO + ^ 


cosh ~ X 
a 


cosh — X 
a 


cos(2/7 + 1) 


^16, 


Deduction 

Case I : If we put Pj = 0 in the eqn. (1 i), we shall obtain Maxwell fluid which is given 
below 


W{x, y,z) = 'Z 

n=0 


(- D" 4P(1 -Xico) 


( 2 / 2 + 1 ) 


nr 


cosh— A' 
a 


cosh — 
a 


e cos( 2/2 + 1) 


%ay 

2b 


where 


K- = [nr ~ a)(l - XiCO)) a" 


Case n : Putting = 0 and = 0 in the eqn. (1 1), we shall obtain purely viscous fluid 
which is given below : 


w(a, y, t) = 


I 


n=0 


i-iT 4P 
(2n+l) Tcr 


cosh— X 
a 


cosh 


K 


e'^cosiln+l) 


nay 

2b 


where = (.n? - (o) 

Numerical Calculation and Discussion 

For numerical calculation of the velocity for Oldroyd first and second order fluid given 
by the eqn. (11) and (16) we take X, = 0.0023, ju, = 0.0005, a = 0.5, b = 0.25, x = 0.75, 
y = 0.45, (0 = 10, m = 10, P = 8, = 0.000023, ti 2 = 0.000005. 

The velocity of a fluid element starts with a maximum value initially and thereafter it 
gradually diminishes asymptotically to a vanishingly small value as time advances and 
tends to infinity. It is represented in the forms of tables and graphs. 
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Fig. 1 - Velocity profile for first order Oldroyd fluid. 



Fig. 2 - Velocity profile for second order Oldroyd fluid 
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Table 1 

Table 2 


Oldroyd fluid 

Oldroyd fluid 


(1st Order) 

(2nd Order) 

t 

W (X, y, t)/A 

W (X, y, t)/A 

0 

24.8646136 

24.9189382 

0.05 

15.0811504 

1 15.1141 

0.1 

9.1471801 

9.167165 

0.15 

5.5480452 

5.5601666 

0.2 

3.3650595 

3.3724115 

0.25 

2.0410117 

2.045471 

0.3 

1.2379362 

1.2406408 

0.35 

0.7508462 

0.7524867 

0.4 

0.4554112 

0.4564062 

0.45 

0.2762209 

0.2768243 

0.5 

0.1675364 

0.1679024 


A General Visco-Elastic Model of Liquid 


A new general empirical model of visco-elastic fluid has been suggested by one of the 


authors (P.R.S.) in the following form 
Xy = - p5^. -t V 


18-20 


\ V > 


V = 2\l 


* I; 9'' 


^ij 




where Xy is the stress tensor, x,y' is the deviatoric stress tensor, the rate of strain tensor, 
p the fluid pressure, Xj are new material constants of which the greatest value Xj represents 
the relaxation time parameter and , A.3 , are additional material constants; \ij are 
also new material constants of which the greatest value pj represents the strain rate 
retardation time parameter and P2 > M'S Mn additional material constants representing 
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the behaviour of a very wide class of visco-elastic liquids the metric tensor in cartesian 
co-ordinates and \i, the co-efficient of viscosity and v,. the velocity components. The material 
constants Xj and jiy designating visco-elasticity satisfy the following conditions 
Xi>X2>X ^> .... > and > P2 > IX3 > .... > p„ i.e. they are arranged in descending order 
of magnitudes. These rheological equations are concerned with the description of the 
mechanical properties of numerous materials under various deformation conditions when 
the materials exhibit the ability to flow the accumulative recoverable deformations. In fact, 
materials satisfying elasticity of solids and viscosity of liquids fall into this category. 
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Abstract 

The object of the paper is to present a theoretical analysis of free convective two dimensional unsteady 
flow with radiation through a porous medium of variable permeabi lity, bounded by an infinite vertical 
porous plate with uniform suction. The permeability of porous medium fluctuates with time about a 
constant mean. Using perturbation method, the expressions for velocity, temperature and skin-friction 
are obtained and the effects of radiation, variable permeability and Prandtl number on velocity, 
temperature and skin-friction amplitude and phase are presented graphically. It is found that the 
velocity decreases with radiation parameter and Prandtl number but increases with permeability. 

Introduction 

Raptis et al} studied the steady free convective flow through porous medium bounded 
by an infinite plate. The same problem, when there is a constant heat flux on the plate, was 
studied by Raptis et al? Thereafter, Raptis^ and Raptis and Predikis*^ studied unsteady flow 
through porous medium, when the temperature of porous plate oscillates in time about a 
constant mean. In above problems, the permeability of porous medium was assumed to be 
constant. In fact, a porous material containing the fluid is a non homogeneous medium and 
the inhomogeneties which can be present in porous medium are numerous. Thus taking 
permeability variations into consideration, Chandrasekhara et al^ and Vedha Nayagam et 
al!" have studied the buoyancy induced flow behaviour adjacent to a horizontal surface in 
a porous medium. Singh et alJ and Singh and Suresh Kumar^ have studied the effects of 
permeability variation on free convective flow through a porous medium bounded by an 
infinite vertical porous plate with constant suction. When the temperature of the plate is 
high, the radiation effects are not negligible. Recently, the effect of radiation and magnetic 
field on the free convective flow along semi-infinite vertical plate bounded by non-porous 
medium has been studied by Takhar et al^ 

The purpose of this paper is to study the effects of radiation and permeability variation 
on free convective flow through porous medium bounded by an infinite porous plate with 
uniform suction. 
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It has been pointed out by Rudraiah et al. that the results of this type of problem 
are of interest in metal production, underground bed flow, oil recovery from partially 
depleted reservoirs, agriculture engineering and various applications in chemical 
engineering. 


Analysis 

Consider the free convective flow of grey gas, near equilibrium in the optically thin 

limit, through a porous medium bounded by an infinite vertical plate with constant suction. 

/ / 

The plate is maintained at constant temperature and is the temperature of the fluid 
far away from the plate. The x' axis is taken along the plate and / axis normal to it. The 
flow in the medium is entirely due to buoyancy force caused by temperature difference 
between the wall and the fluid. Neglecting the viscous dissipation and Darcy’s dissipation 
temrs, the flow can be shown to be governed by the following equations : 








(1) 

du' 
dt' ^ 

du' 

dy' 

-8^(r 

-TJ) 

+ V 

d^u' V , 

0/ k'it') “ 

(2) 

dr ^ , 

dr 

_ j_ 

d^r 

. J- 


(3) 

dt' '' 

dy' = 

' sP ' 

dy'^ ' 

PCp 

, ' dy' 


where u and v are components of velocity along x! and y' directions, g the acceleration 
due to gravity, (3 the coefficient of volume expansion, T ' the temperature. p,v ,k,Cp 
respectively are density, kinemetic viscosity, thermal conductivity and specific heat of the 
fluid at constant pressure, k ' and are permeability and heat flux respectively. 

Following work of other workers^”* the permeability of the porous medium is assumed 
to be of the form 


r(0 = V (1 + e (4) 

where is the mean permeability of the medium, w' is the frequency of fluctuation, t ' is 
the time and e (« 1) is a constant quantity. 
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Cogley et al . ' ' showed that in the flow of grey gas near equilibrium in the opticlly thin 
limit, the radiative heat flux q' is given by 

^ = 4 (r - TJ) I' (5) 


where f = h,. dk , is the absorption coefficient at the wall and is 

•'() di 


Planck 


function. 


The relation (5) has been used by many authors including Tekhar et aO^ and Grief et 


The boundary conditions are 


y' = 0 : z/' = 0, r = TJ 


y' ^ oo : u = 0 ^ T = TJ 


Integration of eqn. (1) gives 


V = - Va 


where Vq > 0 is a constant with dimensions of velocity and negative sign indicates that 
the suction is towards the plate. 

Introducing the following non-dimensional quantities : 

y = y' Vq / V r = f Vq / 4v , w = 4v vv' / Vq , 

,, = ^'/vq, 9 = (r - TJ)I{TJ - TJ). 

the eqn. (2) and (3) in view of (4) and (5), reduce to the following form 

J: ^ ^ ^ r A - - 

4 3/ 3y ^q(i + e e '""^) 


i ^ M = i- ^ 

4 dy Pr 


- F.0 


( 8 ) 
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where 


V p 6 (r ' - r ') 

Gr= \ (Grashoff number) 


and 


Pr = \i.c Ik, 




F = 4vf/p.Cp.v2 


(Prandtl number) 
(Permeability parameter) 
(Radiation parameter) 


The corresponding boundary conditions become 
y = 0:w = 0, 0=1, 
y-^oo:w = 0, 0 = 0. 


(9) 


The partial differential eqn. (7) and (8) are reduced to ordinary one by assuming the 
following series expressions for velocity and temperature. 


w (y, 0 = «o M ^ (y) + . 

9 (y» 0 = 00 Cv) + ^ 01 (y) + 


(10) 

( 11 ) 


Substituting eqn. (10) and (11) in eqn. (7) and (8) and equating the coefficients of like 
powers of G to zero, the following set of ordinary differential equations are obtained : 


Uq + u, 


0 - Y ^0 = - Gr . % 

/Ca 


W| + Wj - 


^iw 1 
^0 




Wi = - Gr . 01 - — 


J_ 

Pr 


%" + %' - F. 00 = 0 


- 


F + 


iw 


0, = 0 


( 12 ) 

(13) 

(14) 

(15) 
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where the primes denote differentiation with respect to y. 

The boundary conditions become 

y = 0 ; Wq = 0- «i = 0, 00 = I, 6i = 0, 

y : Uq = 0, u^ = 0, 00 = 0, 0, = 0 . (16) 

Eqn. (12) to (15) are second order linear differential equations with constant coefficients, 
the solutions of which are straightforward. Solutions of (12) to (15), satisfying the boundary 
conditions (16), are substituted in eqn. (10) and (1 1), to give 


u (y, t) 


Gr 

Pi 





where. 



0 (y, r) = 



a = [Pr + V(Pr2 + 4P/-.P)J/2, 
r, = [1 + V (1 + 4/ko)]/2 , 


(17) 

(18) 

(19) 



It may be noted here that when radiation parameters, F = 0, the above solutions reduce to 
the solutions obtained by Singh and Kumar^. 

Taking the real part of the above solutions, the velocity field can be expressed in terns 
of fluctuating parts as 


w (y, r) = Wq (y) + e {M^ cos wt - M- sin wt) 


( 20 ) 
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+ iMf = Ml Cy) 


M = e-(i + a))-/2 (A-Qsin^ - S cos ^ +5 


= fii*o.)yn (^_c)cos^ - Bsin^ + - Ae'V 


(l+4/*o) + V(l+4/yto)^ + 


■(1+4//^) + V(1 + Alk^ + 


A = 4Gr/ko Pi w, B = 16Gr/ko (16P? + w^) 


C = 4Gr. wIkQ Pj (16P^ + w^) . 


71 

The expression for the transient velocity, for wt = — , is given by 


M (y, 7t/2w) = Ugly) - e M, . 


The fluctuating parts , M,- and the transient velocity for wt = ^ have been plotted 
in Fig. (1 & 2) respectively. 

Skin friction : The skin-friction at the plate in terms of its amplitude and phase can 
now be calculated from eqn. (20), as 


fw f 3m I Gr 


^f~ ^ .,2 


P-Vo Pt 


{a - rj) + e |Al cos {wt -H a) 
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Fig. 2 - Fluctuating parts of velocity profiles. 
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where 


N = N^ + iNi 

^ (A - C) + B (1 + a) ! 2 - Ba , 

N, = Ar^ + ^ - Ca - (A - Q (1+ a) / 2 

and tan a = . 

Results and Discussion 

The velocity profiles for wt = nl2 are shown in Fig. 1. It is observed that near the 
plate (y « 1) the velocity increases rapidly and then decreases slowly away from the plate. 
Also the velocity increases with permeability kQ but decreases with Prandtl number Fr^and 
radiation parameter F. 

It may be noted from Fig. 2 that fluctuating parts of the velocity - remains positive 
and M- remains negative for all values of permeability kQ , Prandtl number Pr, frequency 
w and radiation parameter F. Also, decreases with increase in kQ , Pr, w and F whereas 
M, decreases as Pr or F or w decreases and as kQ increases. 

For the phase and amplitude of the skin-friction, we observe from Fig. 3 that the phase 
tan a remains negative and decreases with increase in permeability but increases with 
Pr and F, This means that there is always a phase lag. Also the amplitude |A^ decreases 
with /cq, Pr and F. 

The temperature distribution is plotted in Fig. 4. The temperature in thermal boundary 
layer decreases with increase in Prandtl number Pr or radiation parameter F. 

Conclusions 


It may be concluded from the above analysis that the effect of radiation is to decrease 
the velocity and temperature in the free convective boundary layer. But, with the increase 
in permeability parameter kQ , the velocity increases in the boundary layer. 
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Abstract 


The photoconductivity of (B^O^ - ZnS) Cu, Cl composites has been studied under various parameters 
such as field, light intensity, temperature and radiation wavelength etc. The (10% BP- ~ 90% ZnS) 
Cu (0-l5%) Cl (1 %) composite shows maximum photoconductivity. The /-V^characteristics are linear 
up to 4V followed by saturation. For moderate temperatures, the photocurrent rises with temperature 
whereas with intensity it saturates. 

(Keywords : photoconductivity/semi-conductors/composites) 


Introduction 

The interaction of photons having energy equal to or greater than the band gap, with 
bound electrons of lattice atoms, creates free electron hole pairs^ The increased carrier 
concentration in presence of light is function of various parameters such as intensity of 
illumination, energy of radiation, applied electric field, temperature^ etc. These effects are 
well understood^ by studying the variation of photocurrent as a function of above parameters. 
A good photosensitive material should not only show a large change in conductivity but 
also respond fast. If trapping centres are in abundance, the response time is slow. Trapping 
also increases the decay time, as the carriers are slowly released after removal of excitation 
source. 

Although a lot of literature exists on the photoconducting properties of single crystals"^'^ 
relatively fewer measurements have been made on binder layers. 

Photoconducing properties of a large number of materials have been investigated by 
several workers^'^. The temperature dependence of photocurrent provides a good 
information about the localized defect states and location of fermi level. 
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The nature of current V 5 . intensity curve gives an idea about the charge trapping and 
recombination processes taking place inside the material. 

In current years properties of gallium-doped ZnO deposited onto glass by spray 
pyrolysis*, optical and electrical characteristics of aluminium doped ZnO thin films prepared 
by Solgel technique^ assembly and applications of the inorganic nanocrystals in polymer 
networks'® etc. have been studied. 

In the present study, an attempt has been made to study the photoconducting properties 
of (B ,03 - ZnS) Cu (0.1%) Cl (1%) binder layers under different parameters to get the 
optimum conditions. 


Materials and Method 

(i) Preparation of sample : For the preparation of sample, the two base materials viz. 
ZnS and B 2 O 3 were taken in different proportions, say. 

(i) 95% ZnS - 5% B 2 O 3 

(ii) 90% ZnS - 10% B 2 O 3 

(iii) 85% ZnS - 15% B 2 O 3 

Then copper impurity of concentration (0.05%, 0.1%, 0.15%, 0.2%) and chlorine (1%) 
was added in the form of NH 4 Cl aqueous solution into the mixture making it a homogeneous 
paste. The mixture was then slowly dried up at about 80“C. It was finely powdered and 
filled into silica tube with graphite corks loosely fitting the ends. For firing, a cylindrical 
furnace, both ends fitted with ruber corks is taken. Temperature was regulated by adjusting 
the external resistance put in the circuit and was measured by a (Cr-Al) temperature 
indicator. The sample was fired at 820'’C for 45 min. under controlled atmospheric 
conditions. 

(ii) Fabrication of cell : Cell was fabricated in the form of parallel plate capacitor by 
using polysterene as binder. The sample suspended in binder was pressed between two 
conducting glass plates. Conducting glass can be made by spraying stanic chloride (SnCl 4 ) 
solution on non-conduting glass plate held at 600T. A thin layer of tin oxide (Sn 02 ) is 
thus deposited which is conducting and transparent. 

(iii) Measurements : For measurement purpose a cell of area -2.5 cm^ with a thickness 
-0.5 mm was mounted in a dark chamber and a d.c. field (50-100 V/cm) was applied and 
light from 300 W Hg lamp was allowed to fall on the surface of the cell. Photocurrent was 
measured through a nA meter under different parameters such as intensity, temperature. 
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field, wavelength of radiation etc. The temperature of cell is changed by changing the 
temperature of oven and intensity by changing the slit width. Optical filters are used for 
changing the wavelength of radiation. 

Results and Discussion 

Photoconducting properties of (5% B 2 O 3 - 95% ZnS), (10% B 2 O 3 - 90% ZnS) and 
(15% B 2 O 3 - 85% ZnS) with different concentrations of Cu and Cl were studied. It was 
found that (10% 8^03 - 90% ZnS) sample with different concentrations of Cu shows 
maximum photoconductivity in comparison to other samples hence the photoconducting 
properties of (10% B 2 O 3 - 90% ZnS) sample with different parameters was studied. 
Photoconductivity changes with the amount of Cu, and is maximum for 0.15% Cu. 
(Fig. 1). 

(i) Ejfect of field : Fig. 2 shows variation of photocurrent with voltage. In presence of 
light initially photocurrent increases with increase in voltage later it saturates. 

The saturation*’^ ‘ of Ip vs. V curves at higher voltage may be explained on the basis 
of class II states, which are imperfection centres lying close to valence band and having 
higher capture cross section for holes than for electrons*“. 

There are also other type of states i.e. class I states which have roughly similar cross 
section for electrons and holes. Under illumination there exist two fermi levels, one for 
electrons and other one for holes. By increasing light intensity the electron fermi levels 
and hole fermi levels are shifted towards conduction band and valence band respectively. 
At a particular voltage and intensity all the class II states may be converted into 
recombination levels. If the voltage is increased beyond these values more electrons will 
be injected into the material. There by raising the fermi level up towards the conduction 
band. The upward movement of the fermi level would convert some of the class II state 
into trapping level causing desensitization with increasing voltage. 

(ii) Effect of temperature : Fig. 3 shows effect of temperature on photoconductivity of 
(10% B 2 O 3 - 90% ZnS) Cu (0.15%) Cl (1%) sample*^. Photocurrent increases with increase 
in temperature-graphs are linear but at higher voltages there is change of slope. 

Fermi levels are shifted towards the middle of the gap with the increase in temperature. 

Arrhenius plots i.e. variation of photocurrent with IIT plots (Fig. 4) are straight lines 
with changing its slope at IIT = 2.9 x lO'^K. Nature of the graphs do not change at 
different supply voltages. 
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Fig. 1 ~ Variation of phtocurrent with % of Cu in the ZnS -* B 2 O 3 composites. 



Fig. 2. - Variation of photocurrent with voltage for 3 composites with 0.15%, 0.1%, 0.05% of Cu. 
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Fig. 3 - Variation of photocurrent with temperature for sample (10% B 2 O 3 - 90% ZnS) Cu (0.15%) Cl (1%). 



Fig. 4 - Arrhenius plots for sample (10% B2O3 - 90% ZnS) Cu (0.15%) Cl (1%). 
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(iii) Effect of intensity : Fig, 5 shows variation of photocarrent with intensity of 
illumination. Initially photocurrent increases with increase in intensity later it becomes 
constant. This means that at a particular intensity all charge carriers are released. 

(iv) Effect of wavelength : Fig. 6 shows variation of photocurrent with wavelength. 
Initially photocurrent decreases with increase in wavelength then becomes constant, The 
dip in phoiocurrent may be due to absorption of radiation by non- radiative centres. 

Conclusions 

The (B2O3 - ZnS) composites shov/ maximum photoconductivity for (10% 6^03 - 
90% ZnS) Cu (.15%) Ci (1%) composition. The 1 - V characteristics are linear upto 4V, 
followed by saturation. For moderate temperature, the photocurrent rises with temperature, 
whereas it saturates with rise in intensity. 
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Abstract 


Water samples from different sources were analysed for physico- chemical and bacteriological 
properties. It was observed that the water available in Gorakhpur City was moderately hard and 
contained appreciable amount of iron. The water available from water works and from tubewell source 
could be safely used for drinking purpose but those from shallow hand pumps and open wells should 
be avoided due to high level of turbidty, iron, chloride and coliform density. 

(Keywords: physico-chemical & bacteriological properties/ water quality.) 

Water is essential for maintenance of all forms of life. It is estimated that water available 
for human consumption is only three percent of total quantity of water available on this 
planet. The demand for water is always increasing as a result of continued population 
growth and rapid agricultural and industrial expansion. It is high time that apart from 
judicious planning and management of this precious resource, the issue of its quality 
determination be seriously taken up for its safe and efficient use. 

The studies on water (Fig. 1) for human consumption have been taken up extensively 
in developed countries but the data for developing countries including India is rather 
scanty^"^^ In the present work, the study on the quality of water for human consumption 
has been taken up in the city of Gorakhpur (Fig. 1) by determining its physico-chemical 
and bacteriological characteristics. The city of Gorakhpour is located in the eastern part of 
State of Uttar Pradesh of India. It is situated on the banks of river Rapti within alluvial 
plain and abundant quality of water is available at shallow depth so that general public 
prefers to draw water from hand pumps and tubewells instead of depending on water 
supplied from city water works which is available in limited localities. 

The quality of water is generally expressed in terms of different physico-chemical and 
bacteriological parameters. The physical parameters of interest are turbidity, colour, taste 
and odour where as pH, alkalinity, chloride, hardness, sodium, potassium, calcium 
manganese, mercury, selenium, iron, arsenic, sulphates, phenolic compounds, nitrates, 
fluoride etc. represent the main chemical characteristics. As far as the question of 
bacteriological characterization goes, coliform density is used as an indicator of faecal 
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Table 1 - 

Details of sample collection. 



Sampling point Location 

Source type 

Sanitary condition 

1. 

Rajendra Nagar 

Deep borewel! 

Good 

2. 

Humayunpur 

Shallow handpump 

Good 

3. 

Fertilizer Colony 

Water works 

Good 

4. 

Medical College 

Mark II handpump 

Bad 

5. 

Shahpur Square 

Shallow handpump 

Fair 

6. 

Dharamshala 

Shallow handpump 

Bad 

7. 

Golghar 

Shallow handpump 

Bad 

8 . 

Head Post Office 

Shallow handpump 

Bad 

9. 

Nakas Chowk 

Mark II handpump 

Good 

10. 

Maya Bazar 

Water works 

Fair 

11. 

Gorakhnath 

Water works 

Fair 

12. 

University Office 

Water works 

Fair 

13. 

Luxhipur 

Open well 

Bad 

14. 

Engineering College 

Shallow hand pump 

Good 

Table 2- 

Details of experimental methodology employed for determining physico-chemical and bacteriological 
parameters of interest. 

S.No. 

Parameter of interest 

Experimental methodology 

1. 

pH 

Glass electrode method 


2. 

Turbidity 

Nephelometric turbidity method 

3. 

Alkalinity 

Indicator method 


4. 

Hardness 

EDTA titrimetric method 


5. 

Chloride 

Agentometric method 


6. 

Sodium 

Flame photometric method 

7. 

Potassium 

Flame photometric method 

8. 

Calcium 

Flame photmetric method 

9. 

Iron 

Phenanthroline method 


10. 

Manganese 

Persulfate method 


11. 

Nitrate 

UV-spectrophotometric method 

12. 

Suplhate 

Turbidimetric method 


13. 

Fluoride 

Alizarin photometric method 

14. 

Coliform density 

Multiple tube fermentation method 
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Table 4 - Drinking Water Standard (See Ref. 16). 

(a) Physical and chemical standards : The physical and chemical water quality of water should not 
exceed the limits shown in the table below: 


S.No. 

Characteristics 

Acceptable* 

Cause for Rejection** 

1 

2 

3 


4 


1. 

Turbidity (Unit on NTU Scale) 

2.5 

(5) 

10 

(10) 

2. 

Color (Unit on Pt-Co Scale) 

5.0 


25 


3. 

Taste and odor 

unobjectionable 

Unobjectionable 

4. 

pH 

7.5 to 8.5 

(6.5-8.5) 

6.5 to 9.2 (NR) 

5. 

Total dissolved solids(mg/l) 

500 

(400) 

1500 

(2000) 

6. 

Total hardness (mg/1) as 

200 

(300) 

600 

(600) 

7. 

Chloride (as Cl) mg/1 

200 

(250) 

1000 

(1000) 

8. 

Sulfate (as SO 12) 

200 

(200) 

400 

(400) 

9. 

Fluorides (as F) mg/I 

1.0 

(1.0) 

1.5 

(1.5) 

10. 

Nitrates (as NO3) mg/1 

45 

(45) 

45 

(100) 

11. 

Calcium (as Ca) mg/l 

75 

(75) 

200 

(200) 

12. 

Magnesium (as Mg) mg/I 

30 

(30) 

150 

(100) 

13. 

Iron (as Fe) mg/l 

0.1 

(0.03) 

I.O 

(1.0) 

14. 

Manganesf' (as Mn) mg/l 

0.05 

(0.03) 

0.5 

(0.3) 

15. 

Copper (as Cu) mg/l 

0.05 

(0.05) 

1.5 

(1.5) 

16. 

Zinc (as Zn) mg/l 

5 

(5) 

15 

(15) 

17. 

Phenolic compounds 

0.001 

(0.001) 

0.002 

(0.0032) 


(as Phenol) mg/l 





18. 

Anionic detergents (mg/l) 

0.2 

(0.2) 

1.0 

(1.0) 

19. 

Mineral oil (mg/l) 

0.01 

(0.01) 

0.3 

(0.3) 

Toxic Materials 





20. 

Aresenic (as As) mg/l 

0.05 

(0.05) 

0.05 

(NR) 

21. 

Cadmium (as Cd) mg/l 

0.01 

(0.01) 

0.01 

(NR) 

22. 

Chromium (as Cr’ j mg/I 

0.05 

(0.05) 

0.05 

(NR) 

23. 

Cyanide (as CN) mg/l 

0.05 

(0.05) 

0.05 

(NR) 

24. 

Lead (as Pb) mg/l 

0.1 

(0.05) 

0.1 

(NR) 

25. 

Selenium (as Se) mg/l 

0.01 

(0.01) 

0.01 

(NR) 

26. 

Mercury (total as Hg) mg/l 

0.001 

(0.001) 

0.001 

(NR) 


Values in parenthisis are the Ref. 16. 

Values IS- 10500, 1991 with January, 1993 amendments. 

(b) Bacteriological standards : Coliform count in any sample of 100 ml should be zero. 

Water in the distribution system shall satisfy all the three criteria indicated below : 

(1) £. Coli count in 100 ml of any sample should be zero. 

(2) Conforms organism not more than 10 per 100 ml shall be present in any sample. 

(3) Coliform organism should not be more detectable in 100 ml of any two consecutive samples or more than 
50% of the samples collected for the year. 

*The figure indicated under the column acceptable are the limits up to which the water is generally 
acceptable to the consumers. 

^Figures in excess of those mentioned under ‘acceptable’ render the water not acceptable but still may be 
tolerated in the absence of alternative and better source but upto the limits indicated under column “cause for 
rejection” above which the supply will have to be rejected. 
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pollution rather the pathogens themselves. The assessment of microbial safety of water, 
therefore, is expressed in terms of coliform (MPN) in 100 ml of water. 

A total set of five samples from fourteen locations selected at random were collected 
and analysed for various water quality parameters. Table 1 gives detail of the sample 
collected in terms of the location of the source, the type of source and the sanitary conditions 
prevailing around the source. If the area near the source was clean and hygenic, the sanitary 
condition was termed good and in case a platform was provided near the source but drainage 
facility was not adequate, the sanitary condition was termd fair. The term bad condition 
was given to those sources which were not provided with platform and had very bad 
drainage facility or a big public drainage was passing through nearby area. 

The water samples were collected in plastic bottles. The bottles were first rinsed with 
distilled water and then with sample water before actual collection. Water from the source 
of interest was allowed to run for two min. before taking samples to ensure flushing of 
stagnant water if any in the domestic water supply pipes or the handpumps. For open well, 
the samples was collected by dipping the sampling bottle tied with string taking care all 
along that the bottle did not touch the sides of the well. The water samples for bacteriological 
analysis were collected in 100 ml sterilised glass bottles. 

The physico-chemical and bacteriological analysis was carried out by employing 
standard methods^"^ using reagent which were prepared and stored exercising all precautions 
mentioned in the above stated standard methods. Table 2 gives details of the experimental 
methodology used for determining the water qualiy parameters of interest. 

The arithmetic mean values and the associated standard deviation for the 
physico-chemical and bacteriological parameters of interest are shown in Table 3. The 
values were obtained from five sets of samples which were taken in the period between 
22nd October 1996 and 8th January 1997 at intervals of two weeks. (Table 4). 

The pH values of water samples analysed in the present work varied from 7.6 to 8.3 
and were found to be within permissible limit. Turbidity which reflects on the appearance 
of water may not be important from health point of view but plays an important role in 
consumer satiafaction. The turbidity of typical treated water is expected to range between 
0 to 2.5 NTU. In the present study, the tutbidities of water samples analysed varied from 
1.2 to 7.2 NTU, the highest value was obtained from water sample taken from an open 
well and lowest from water samples taken from deep bore well. 

Alkalinity which is the capacity of water to neutralize strong acid, itself is not harmful 
for human beings but the desired value for domestic water supply is expected^^ to be less 
than 100 mg/1^^. In the present work, alkalinity of water samples ranged from 104.2 to 
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271.6 mg/1, the highest value being shown from water samples taken from shallow hand 
pump source. Hardness which understood in terms of capacity of water to precipitate soap, 
was found to vary between 67.6 and 150 mg/1 of equivalent calcium carbonate thereby 
indicating that the water available in the city of Gorakhpur could be categorised as 

moderately hard. 

The chloride content which generally reflects the presence of sodium chloride and is 
some times used as an indicator of pollution give a noticeable salty taste to water if its 
concentration rises above 250 mg/I. In the present work, chloride concentration ranged 
from 25.2 to 460 mg/1. The high values were obtained again from water samples taken 
from shallow hand pumps. 

The concentration of sodium, potassium, calcium, manganese, sulphates, nitrates and 
fluorides were also tested and were observed to be within the permissible iimits^^’^^. 
However, the concentration of iron was observed to be on the higher side. Iron though 
present in biologically significant proteins, is not a cause for concern as far as the general 
health of population goes but it has a nuisance value since it precipitates as ferric hydroxide, 
a rusty coloured silt which makes water unpalatable and causes stain to the laundry and 
plumbing fixtures. 

The most common and widespread danger associated with drinking water is 
contamination either by sewage, other waste or human and animal excrement which have 
been known to cause water borne diseases. The universally acceptable principle for 
monitoring and assessing the microbial safety of water supply is through the measurement 
of conform density (MPN). In the present work, the water samples from open well and 
shallow hand pump showed high level of coliform density. 

In the present work a modest attempt was made to evaluate the standard of drinking 
water available to the citizens of Gorakhpur City and it was observed that the water available 
to them was moderately hard and had appreciable concentration of iron. It is recommended 
that people should draw water from sources having deep boring but should avoid drawing 
of water from open wells and shallow hand pumps. 
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